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We study via RG, numerics, exact bounds and qualitative arguments the equilibrium Gibbs 
measure of a particle in a d- dimensional gaussian random potential with translationally invariant 
logarithmic spatial correlations. We show that for any d > 1 it exhibits a transition at T = T c > 0. 
The low temperature glass phase has a non trivial structure, being dominated by a few distant states 
(with replica symmetry breaking phenomenology). In finite dimension this transition exists only in 
this "marginal glass" case (energy fluctuation exponent 9 — 0) and disappears if correlations grow 
faster (single ground state dominance 9 > 0) or slower (high temperature phase). The associated 
extremal statistics problem for correlated energy landscapes exhibits universal features which we 
describe using a non linear (KPP) RG equation. These include the tails of the distribution of 
the minimal energy (or free energy) and the finite size corrections which are universal. The glass 
transition is closely related to Derrida's random energy models. In d = 2, the connexion between 
this problem and Liouville and sinh-Gordon models is discussed. The glass transition of the particle 
exhibits interesting similarities with the weak to strong coupling transition in Liouville (c = 1 barrier) 
and with a transition that we conjecture for the sinh-Gordon model, with correspondence in some 
exact results and RG analysis. Glassy freezing of the particle is associated with the generation 
under RG of new local operators and of non-smooth configurations in Liouville. Applications to 
Dirac fermions in random magnetic fields at criticality reveals a peculiar "quasi-localized" regime 
(corresponding to the glass phase for the particle) where eigenfunctions are concentrated over a 
finite number of distant regions, and allows to recover the multifractal spectrum in the delocalized 
regime. 



I. INTRODUCTION 

Despite significant progress in the last two decades, dis- 
ordered systems continue to pose considerable theoretical 
challenges. Two important questions still largely open, 
are, respectively, to which extent the (better understood) 
mean field models are relevant to describe low dimen- 
sional physical systems, and, in the special case of two 
dimension, to which extent the powerful field theoretic 
treatments developed for pure models can be adapted to 
treat disordered models. 

A celebrated controversy is whether the structure 
found in the solution of mean field models for spin glasses 
and other complex disordered systems, both in the statics 
jlj and in the dynamics has any counterpart in the 
world of experimentally relevant low dimensional mod- 
els. Specifically it has been vigorously questionned || 
whether the breaking of the phase space in "many pure 
states", predicted to occur in mean field, may also oc- 
cur in short range models, and how it can be properly 
defined ^J^j . The unusual nature of the technique used 
to solve the statics, i.e the replica method with a hier- 
archical breaking of the permutation symmetry between 
replicas in the limit n — > (RSB), did not contribute 
to make the physics transparent. A distinct structure, 
which remarkably parallels the one in the statics, has 
been found Q to occur in the nonequilibrium dynamics. 
The dynamical problem can be studied by a priori bet- 
ter defined methods and leads to predictions which are 
in principle directly testable in experiments, such as a 



non trivial generalization of the fluctuation dissipation 
relations. Even so, it has been emphasized that mean 
field models, which usually involve infinite range or infi- 
nite number of component limits, may not capture phys- 
ical processes important in low dimensions. The alterna- 
tive "droplet picture" in its simplest form || postulates 
the existence of a single ground state with excitations 
(droplets) of (free) energy AE scaling with their size x 
as AE ~ x e , 9 > 0. It provides a more conventional scal- 
ing description of the glass physics, as being controlled 
by zero temperature RG fixed points where temperature 
is formally irrelevant (with eigenvalue —0). 

Another important advance was the exact solution of 
simpler prototype models, such as the random energy 
model (REM) || , where one consider simply a collection 
of independently distributed energy levels, as well as its 
generalization, the GREM 0, or the Directed Polymer 
on the Cayley Tree (DPCT) with disorder g. These so- 
lutions being direct with no use of replica, their results 
can be fully relied upon. They exhibit a similar physics 
though, with a glass transition and in the glass phase, 
an exponential tail for the distribution of the free energy 
P(f) ~ e>3c * f° r negative /. This feature is crucial to 
recover the same physics, and indeed many observables 
were found to be similar j^] . In fact the alternative solu- 
tion of the REM using replicas, given in § or the one of 
the DPCT do involve RSB. In the REM model the 
structure of the glass phase is particularly transparent as 
being dominated by a few states H|| . 

It is important to go beyond models defined in mean 
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field or on hierarchical (or ultrametric) structures and to 
study simple yet non trivial (and non artificial) finite d 
models with full statistical translational invariance. In 
this paper we study the model of a particle in a gaussian 
random potential V(r) with spatial correlations which 
are invariant by translation and which grow as the log- 
arithm of the distance. We consider this model in any 
dimension d, but in d = 2 it has also been studied re- 
cently since it is of direct relevance for several physical 
systems [|ll|-]l9||. One example is a spin model with XY 
symmetry and random gauge quenched disorder, which 
arises naturally in describing Josephson junction arrays 
p0| or 2D crystalline structures with smooth disorder, 
e.g. flux lattices in superconductors pi] , or electrons at 
the surface of Helium p^| . In this model, a single topo- 
logical defect (a XY vortex), or a single neutral pair, 
sees precisely a random potential with logarithmic cor- 
relations J14|-p8[. Another example arises in a model of 
localization of Dirac fcrmions in a random magnetic field, 
motivated by quantum Hall physics. There, the zero en- 
ergy E = normalized wavefunction is identical to the 
Boltzmann weight of the particle studied here [pl]-p^|. 
This wave function is "critical" in a sense discussed be- 
low. 

Here we study this model using a renormalization 
group (RG) approach, bounds, numerical methods and 
qualitative arguments. We show that it exhibits a tran- 
sition at T = T c > in any d > 1. We find that in the 
high temperature phase the particle is essentially delocal- 
ized over the whole system, while in the low temperature 
glass phase the Gibbs measure is concentrated in a few 
minima. The fact that such a simple (finite d) model 
exhibits a genuine glass transition is already surprising. 
Indeed, as we argue, this transition exists only for such 
a "marginal" type of correlations (which correspond to 
9 = in the glass scaling mentionned above |23|]). It 
disappears (for gaussian V(r)) if correlations grow faster 
(with only a low temperature phase and single ground 
state dominance) or slower (with only a high temperature 
phase) . Logarithmic growth of correlations thus produces 
exactly the right balance between the depth of the energy 
wells and their number (entropy). Note that for slower 
growing correlations one can recover a transition but only 
by artificially rescaling the disorder variance with the size 
of the system: in the extreme case of uncorrelated vari- 
ables, it is the REM model. Here by contrast, there is 
a genuine phase transition in the thermodynamic limit, 
with no need for rescaling. Most interestingly, the glass 
phase is non trivial. The Gibbs measure is concentrated 
in a few distant minima which remain in a finite number 
in the thermodynamic limit. This is because the extrema 
of random variables with such correlations exhibit an in- 
teresting property of "return near the minimum": there 
is, with a finite probability in a sample of size L — > +oo, 
at least one second minimum far away (at distances of 
order L), and with a finite energy difference with the 
absolute minimum. And there are not too many (a ther- 
modynamic number) of these secondary minima, leading 



to a zero entropy. As in the REM, this property leads 
here naturally to a non trivial ground state structure, 
reminiscent as we discuss of a genuine property of replica 
symmetry breaking in the replicated theory. The low 
temperature limit corresponds to a non trivial problem 
of extremal statistics of correlated variables, studied here. 

Another interesting property of this model is its re- 
lation to the Liouville model (LM) and the sinhGordon 
model (SGM) in d = 2 (and their boundary restriction in 
d = 1): V(r) turns out to be the Liouville field while the 
LM and SGM partition functions arise simply as gen- 
crating functions of the probability distribution of the 
partition sum Z[V] = f e~P v ^ of a single particle. The 
high temperature phase for the particle corresponds to 
the weak coupling regime for the LM and SGM, where 
we find that known exact results compare well with re- 
sults for the particle. In the SGM we predict here the 
existence of a transition (more appropriately, a "change 
of behaviour" ) . It corresponds to the glass transition for 
the particulc, which also exhibits interesting similarities 
with the weak to strong coupling transition in the Li- 
ouville theory (and the so called c = 1 barrier). The 
glassy freezing of the particule is associated, in the LM 
and SGM, to new local operators and non smooth con- 
figurations being generated under RG. 

To study the model we will introduce a RG approach 
based on a Coulomb gas renormalization a la Kosterlitz. 
It leads to a non linear RG equation (of the so-called 
Kolmogorov KPP type) for the full probability distribu- 
tion of the "local disorder" . Indeed, a separation be- 
tween the long range part of the disorder and the local, 
short range part, arises naturally in our approach. The 
RG equation has traveling wave types of solutions. The 
corresponding well known problem, in such non linear 
equations, of the selection of the velocity of the traveling 
front, and its freezing for T < T c is related to glassy freez- 
ing of the particle free energy and, in the LM or SGM, 
to the "selection" of the anomalous dimensions (and at 
the transition dimension degeneracy leads to logarithmic 
operators) . When temperature is lowered the local disor- 
der becomes broadly distributed and the freezing occurs 
when its tails become relevant. Our RG method indi- 
cates that the physics depends only weakly of d. We will 
take advantage of this fact and check our results using 
simulations in d = 1. 

It is important to compare the present work to pre- 
vious studies of the model. The existence of a freez- 
ing transition in d = 2 has been conjectured previously 
0,[5|Jl6) based on an approximation which completely 
neglect spatial correlations (REM approximation, see be- 
low) . Stronger arguments were given in |l3| , but did not 
fully establish the existence of a transition, which is done 
here (see Appendix @). The present work is thus mo- 
tivated by the need to go beyond the REM approxima- 
tion to describe this problem. In particular one wants to 
know what is the precise universality class of the model, 
which we hope can be determined from the RG method 
introduced here. This RG method yields some remark- 



2 



able universal features of the probability distribution of 
the free energy and of its finite size corrections, differ- 
ent from the REM. It shows that the problem is more 
closely related to the directed polymer on a Cayley tree. 
A qualitative analogy between the present model and 
the DPCT was in fact cleverly guessed recently in Ref. 
jl6],[ll]]. It is based on the observation that the energy 
of polymer configuration on a tree also scale logarithmi- 
cally with the overlap distance defined on the tree (see 
Fig [l]) . It is remarkable that this connection naturally 
emerges here from the Kostcrlitz type RG performed on 
this problem, via the KPP equation. It is all the more 
surprising, since the model studied here has statistical 
translational invariance, while a tree has a hierarchical 
structure. The solution of Derrida and Spohn || (and 
the mapping onto the DPCT proposed in jlffl,nT|) would 
be exact for random variables V^r) correlated with a hi- 
erarchical (i.e ultrametric) matrix of correlation. Here 
instead the correlations are translationally invariant and 
it is thus important to understand the origin of the anal- 
ogy with the DPCT and to which extent it holds. The 
RG procedure developed in this paper is a first attempt 
to address these questions. The result is that we can 
make the mapping precise: at least for the universal ob- 
servables studied here (e.g. the tails of the free energy 
distribution) , the mapping is onto a continuum branching 
process, i.e a continuum limit of a Cayley tree (whereas 
^6| , p| could not be so specific). 




FIG. 1. Directed polymer approximation: sites are the tree 
endpoints. The are uncorrelated of variance 2a. The ran- 
dom potential at i is Vi 



Vj = + v%> + v^>. Thus (Vi - Vj 



Dj 1 ^ + v[ 2 ' + Uj 3 ' and at j it is 



(3) 



4ad(i,j) where 



d(i,j) ~ In \i — j\ is the distance (in generations) on the tree 

The present model has also been studied in the con- 
text of random Dirac problems and localization. An 
early study [M of the E — wavefunction established 
that it was critical (in the sense of corresponding to a 
"delocalized" wave function, while E ^ has finite lo- 
calization length). However this study missed the glass 
transition. Later studies |]l3f computed the multifractal 
spectrum based on the REM approximation and noticed 
the existence of a strong disorder regime. These and 
other studies [ pdp^ ] however focused on properties of 
the high temperature phase: it was conjectured |l2j that 
the (conformal) Liouville field theory (LFT) (i.e a con- 
tinuum limit of the LM) was able to describe all spatial 



correlations of the model in the high temperature phase. 
These works call for more investigations. First the glass 
transition and the peculiar physical properties of the low 
T (i.e strong disorder) phase have not been addressed, 
even at the most qualitative level. We thus find it use- 
ful to present the problem from a different perspective 
by comparing with other types of correlated disorder, or 
by recasting it as a problem of extremal statistics. Al- 
though well known properties [^5| of extremal statistics 
of uncorrelated variables were often used to study model 
disordered systems (see e.g. fp6|-p8[) a lot remains to be 
understood about the (more realistic) case of correlated 
variables. Second, the question of the universality class is 
in our opinion far from established. Evidence for the LFT 
description mostly comes from reproducing the multifrac- 
tal spectrum as given by the REM approximation and one 
would like to check it against more detailed predictions. 
The present RG procedure is a step towards clarifying the 
connexion between this model and solvable models such 
as Derrida's REM and Derrida-Spohn DPCT. In this re- 
spect finite size corrections are important to understand, 
as they are found to exhibit universal prefactors allow- 
ing to distinguish between various universality classes. 
In addition, they determine the anomalous dimensions, 
and thus control the critical behaviour, in the full disor- 
dered XY model as shown in p8[ . Since they are found 
to be very large, they are also crucial in order to ana- 
lyze the results of numerical simulations. In particular, 
although we confirm the result of we also conclude 
that the sizes used in the numerical study of fll| were in 
fact vastly insufficient for drawing firm conclusions: we 
do perform here a more detailed finite size analysis on 
much larger samples to confirm analytical predictions. 

The model studied here is thus related to a surprising 
number of interesting problems. Let us mention for com- 
pleteness that it also has connexions to problems such 
as two dimensional interfaces, or films, confined between 
two walls (for [3 = +oo it is the confinement entropy 
of a film) , wetting transitions |^9| , extremal statistics of 
correlated variables useful e.g. for problems of "persis- 
tence" in nonequilibrium dynamics |3C|| , and finally, to 
the clumping transition of a self gravitating planar gas 
fnfl . We will not explore these connections in details 
here. 

This paper is organized as follows: in sectio n [Ti|, the 
single particle model is defined and in section [I B the 



random energies approximation (REM) is applied, which 
amounts to neglect the spatial correlations of the ran- 
dom potential. The full problem, with correlations taken 
into accou nt, is related to the description of extremal 
statistics in [II C| , and three different classes of correlations 
are identified in [[I D| from qualitative arguments. A new 
renormalization (RG) technique is applied to this prob- 
lem in section III. The resulting non-linear scaling equa- 
tion for the d istribution of the local disorder is studied in 
section III C , and found to be related to the Kolmogorov 
KPP equation, which admits front solutions. This con- 
nection between front solutions of non-linear equations 
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and the r enorm alization of disordered models is pursued 
where a solution to the REM is found via 



m section 



hid 



a similar non-linear RG (details in Appendix C). The non 



trivial nature of the glass phase is discussed in III E to- 
gether with its relations to replica symetry breaking. In 
part IV we present a numerical analysis of the problem 
of the particle in a random potential in d = 1. Section ^ 
is devoted to the connexion between the particle model 
-and its transition- and entropic phenomena in the Liou- 
ville and Sinh-Gordon models. A direct RG analysis in 



Section VC allows to recover the corresponding change 
of behaviours in these models. Section VI contains the 
applications to the properties of the critical wave func- 
tion of a Dirac fermion in a random magnetic field, in 
particular the multifractal properties and the property 
of quasi-localization. Appendix A contains an outline of 
the proof of the existence of a transition, Appendix B 
is a review of well known (and not so well known) re- 
sults about extremal statistics, Appendix D contains an 
extended model which exhibits three phases. 



II. MODEL AND QUALITATIVE ANALYSIS 

In this Section we define the model of a single parti- 
cle in a correlated random potential. Then we describe 
the random energy model (REM) approximation used in 
previous studies which consist in neglecting correlations. 
We then pose the new questions which we want to ad- 
dress here for the true model and present a qualitative 
analysis showing physically why we expect that logarith- 
mic correlations (as opposed to faster growing or slower 
growing correlations) is the only case which leads to (i) a 
glass transition (ii) a low temperature phase with a non 
trivial structure of quasi degenerate distant minima 



A. the model 

The equilibrium problem of a single particle in a d- 
dimensional random potential is defined by the canonical 
partition function 



-f3V(r) 



(1) 



where /3 = 1/T is the inverse temperature, in a sample 
of finite size L (and total number of sites L d ) and for a 
given configuration of the random variables V(r). The 
equilibrium Gibbs measure, or probability distribution 
for the position of the particle is: 



p(r) 



-PV(r) 



/Z[V] 



(2) 



We are interested here in cases where the random vari- 
ables V(r) can be correlated. As discussed below, the 
statics (and dynamics) of this problem in the limit of 



large sizes depends on the type of correlations, the dis- 
tribution of the disorder and the dimensionality of space 
d. Some of these cases and their dynamical aspects (such 
as the Sinai model) have been extensively studied, e.g. 
in the context of diffusion in random media J32 ]. Even 
logarithmic correlations in d = 2 were studied then [j33| , 
but it was not realized at that time that a static glass 
transition could exist in that case. 

Correlated random potentials V(r) are most conve- 
niently studied for Gaussian distributions, on which 
we focus, parametrized by the correlator T(r, r') = 
V(r)V(r') (and we choose V(r) = 0). Non gaussian ex- 
tensions will be mentionned. Unless specified otherwise 
the correlations will be chosen translationally invariant 
r(r,r') = r^(r — r') with cyclic boundary conditions, 
or in (discrete) Fourier space V(q)V(p) = T(q, p) = 
r(q)<5 p ,_ q . We will often denote (V(r) - V(r')) 2 = 
f(r - r') = 2/ q r(q)(l - cos(q.(r - r'))) (with / q = 
1 — > f 

L l^q J (2TT) d >- 



One important quantity is the free energy: 
F[V] = -T\nZ[V] 



(3) 



and, since it fluctuates from configuration to configura- 
tion, as F[V] = F[V] + 6F[V] we will be interested in 
its average F = F[V] and in its distribution. From the 
convexity of the logarithm follows the well known exact 
bound for F in terms of the annealed free energy Fa'- 



-TlnZ = F>F A = -T\nZ 
1 



F>-(TdlnL + — U(r) 2 ) 



(4) 
(5) 



for the Gaussian case. 

In this paper we will mainly focus on the case of cor- 
relations growing logarithmically with distance: 



(V(r) - U(r')) 2 ~ 4a In 



r — r 



a< |r-r'| < L (6) 



which also requires a small distance ultraviolet (UV) cut- 
off a (we can set here a = 1 in accordance with the defi- 
nition |] of a discrete model, but in the following Sections 
we will consider a continuum version and vary a). This 
behaviour is achieved in d dimension by choosing a prop- 
agator in Fourier space T(q) ~ 2 g^ 2 ^ . The d = 2 case 
is also of special interest as the propagator is the usual 
Coulomb one: 



r(q) 



Ana 

„2 



(7) 



and boundary conditions must be specified later on. It is 
importa nt to no te that for LR correlations the single site 
variance U(r) 2 = Ti(0) diverges with the system size, 
e.g. for (§) one has T L (0) - 2abx(L/a). 

For such logarithmic correlations (as well as for weaker 
correlations |54|) one will find that F scales as d In L (con- 
sistent with the number of degree of freedom being L d in 
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this problem). Thus it is natural to define the intensive 
free energy: 



f(P) 



F[V] 



l^+oo d\nL 



lim 



(8) 



which will be found to be self averaging. The above 
bound gives: 



(9) 



Thus we will find that F[V] ~ f{[3)d\nL with subdomi- 
nant corrections. These corrections have a non fluctuat- 
ing universal 0(ln(lnL)) piece, as well as an O(l) fluc- 
tuating part [V] which we will both study. 



B. the REM approximation 

A useful approximation to the problem studied here, 
which can be called the REM approximation, consists in 
neglecting all correlations but keeping the on site vari- 
ance exact |fL7j[L3| : 



r £ (r) 



iREM 



(r) =T L (0)5 r . T/ =2(7 In [ - 
1 a 



(10) 



The corresponding Gaussian REM model can then be 
solved, being identical to and one finds that it ex- 
hibits a transition at (3 C = \fdfa with: 



fW) 



m = ~K 



f3> Pc 



P<(3 C 



(11a) 
(lib) 



Most previous studies of the original model (all in d = 2) 
amount to study the REM approximation and argue that 
it is a good approximation. Indeed, as we will also find 
here, this REM approximation appears to give the exact 
result for some observables (e.g. for /(/?)). In particular, 
it does seem to give correctly the transition temperature 
/3c 



C. beyond the REM approximation: extremal 
statistics of correlated variables 

Since it is not obvious a priori why logarithmic cor- 
relations can be considered so weak as to be neglected, 
one would like to go beyond the REM approximation and 
describe the effect of the neglected correlations . One 
would like to understand why this approximation works 
for some observables (and for which ones) and whether 
it gives exactly the universality class of the model (i.e all 
universal behaviour of observables). The answer to the 
latter is negative: as our analysis will reveal, the correla- 
tions do matter for the more detailed behaviour and the 



original model (|l]J^) is not in the same universality class 
as the REM model. 

In fact, the problem at hand is related to describing 
universal features of the extremal value statistics for a 
set of correlated random variables. Indeed, the zero tem- 
perature limit (T = for fixed L) of the problem defined 
by (Q) amounts to finding the distribution of the mini- 
mum — limy—^o Tin Zl — V m i n = min r ({V r }) of a set of 
correlated random variables. In the case of uncorrelated 
(or short range correlated) variables a lot is known in 
probability theory on this problem (see e.g. pfifl ), some 
of which is summarized in Appendix [b|. For the type of 
distributions considered here (gaussian and some exten- 
sions) the distribution of the minimum V m in has a strong 
universality property, being given, up to - non trivial - 
rescaling and shift (see Appendix [b| and below), by the 
Gumbell distribution: 



Proba(y < x) = Q(x) = exp(— e x ) 



(12) 



The Gumbell distribution thus appears as the distribu- 
tion of the zero temperature free energy in the REM. 
For the case of a Gaussian distribution the standard 
probability theory results are usually given in terms of 
a variable X r such that X% = 1. One can simply rescale 
V r — V2(T In LX r from Appendix [B^ and get: 



V m i n = 2\Tod\nL — ln(47rdlnL) 




(13) 



where y is distributed as in (|l2]). 

Much less is known in the case of variables with 
stronger correlations studied here, though it is more im- 
portant in practice. The statistics of V m in m the loga- 
rithmically correlated case is thus one of the open issues 
discussed here. One key question is to determine what is 
universal in the distribution of the minimum of correlated 
variables. Here, we can formulate the question as follows: 
given Gaussian random variables satisfying @, what in 
the distribution of the minimum (i.e of the ground state 
energy for fixed large L) is universal, i.e depends only on 
a and not on the details of the correlator F(r) at short 
scale. Writing 



Vrr. 



sv m 



(14) 



one finds, for the logarithmic correlator, that the aver- 
aged ground state energy must satisfy: 



Vmiri > —2y adhxL 



(15) 



which follows from the above annealed bound, together 
with the fact that df/dT = —S < 0. Furthermore one 
will find here that Vmin ~ e m i n d In L up to a positive sub- 
dominant - universal - piece and that e m i n = —2y/o~/d 
saturates the bound. In the distribution of 5V m i n ~ 0(1) 
we can clearly expect less universality that in the problem 
of random variables with short range correlations fl37|] . 
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D. Qualitative study of a particle in a random 
potential 

Before describing the RG method which allows to go 
beyond the REM approximation, let us give some sim- 
ple qualitative arguments and numerical results which 
illustrate the main physics of the thermodynamics of a 
particle in a correlated random potential. To put things 
in context we discuss several types of correlations (short 
range, long range, and marginal). We focus on d = 1 for 
simplicity but the arguments extend to any finite d. 

Whether there is a single phase or not here comes sim- 
ply from whether the entropy of typical sites wins or not 
over the energy of the low energy sites. When there is a 
low temperature phase, to decide its structure one must 
pay special attention to distant secondary local minima. 

Indeed, when there is a low temperature phase, it is 
controlled by the regions with most negative potential. 
To investigate its structure one can start, for a given sys- 
tem of size L, with the T = state which is determined 
by the absolute minimum over the system, denoted V m in 
and located at r m i n . At T very small but strictly posi- 
tive, each (low lying) secondary local minima V will also 
be occupied with a probability ~ e -(V-v min )/T w j 1 j c ] 1 j s 
very small except when V — Vmin ~ 0(T). Thus to char- 
acterize the low temperature phase we need to know how 
many of these secondary minima exist and where they 
are located. For a smooth enough disorder (see e.g. Fig 
^) there will always be "trivial" secondary local minima 
in the vicinity of r m i„. To eliminate these, we define 
V m i n 2(R) as the next lowest minimum constrained to 
be at a distance at least R of the absolute minimum. 
An interesting quantity to study is then the distribution 
P R , L (AE) of AE(R, L) = V min2 {R)~V min over environ- 
ments (which a priori depends on R and L) 



P(AE) 



equivalently Tl(t) — > at large r (or e.g. T q ~ q 



-d+S 






P(AE) 




P(AE) 




AE 




: ae 




^ AE 



Q?LJ 



1 



FIG. 2. Three cases for the distribution of energy difference 
AE between absolute and secondary minimum (separated at 
least by R ~ L c ) in a system of size L: (a) short range cor- 
related potentials, AE typ — > logarithmically with size (b) 
algebraically growing correlations AE typ — * +oo (c) logarith- 
mic correlations, AE typ remains constant as the system size 
increases. 

We now distinguish three main cases, according to the 
behaviour of the correlator (V(r) — V(r')) 2 = T(r — r') 
at large scale (we restrict to gaussian potentials [Q). 
In these three cases the distribution Prl(AE) has 
markedly different behaviours as illustrated in Fig. ^: 

(i) short range correlations, i.e T(r) — > Cst at large r, 



with S > 0). In this case it is clear that there is only a 
high temperature phase in any finite d and no phase tran- 
sition. The entropy TdlnL of typical sites (of energy 
typically ~ 0(1)) always wins over the energy of opti- 
mal sites (Vmin ~ \j2o~d\n.L for gaussian distributions 
with on-site variance a ). The optimal energy V m in can 
be estimated using \/L d = P 1 (V)dV in terms of 

the single site distribution Pi(V), which yields the ex- 
act leading behavior for uncorrelated disorder [^9| (and 
also for weak enough correlations - see Appendix |^). 
Thus, the particle is delocalized over the system for all 
T > 0. One estimates the number of states within AE 
of the minimum as N(AE) ~ L d f^+ AE p l {y)dV ~ 

exp(A_E\/2<iln L / \fcr) for a Gaussian distribution. Thus 
there is a large number of sites almost degenerate with 
the absolute minimum V m i n , separated by finite barriers, 
and AE typ decays to as a power of 1/lnL (1/VhiL for 
a Gaussian) EO ]. These minima however are irrelevant 
for the thermodynamics of the system at a fixed finite 
temperature. 

For these minima to play a role and to obtain a tran- 
sition even for SR disorder one needs to perform some 
artificial rescaling, as in the REM model |{|: either, 
at fixed size, to concentrate on the very low T region, 
(e.g. take (3 ~ InL in the Gaussian case), or equiva- 
lently, to rescale disorder with the system size. By mak- 
ing disorder larger as the system increases, for instance 
using Pi(y) - e -\VIVtvv\ a with Vum ~ (lni) 1_1 / Q one 
recovers artificially a transition P8fl . For a — 2 and 
uncorrelated V(r) this is exactly the REM studied in 
||. There, the simple argument for the transition is 
that the averaged density of sites at energy E = V is 
-E 2 /(2a L ) j^fi^l (related to the annealed 
,e-^0(£)). 



L" 



n(E) 

partition sum via Z = J E e~ p ^Q(E)). If = a is not 
rescaled the average energy is O(l) and the huge entropy 
of these states always wins. If a scales with L as ol ~ 
2<T In L then there is a transition at C = \J dja. Indeed, 
n(E) ~ exp(cilni(l - (e/e mm ) 2 )) where e = E/{d\nL) 
and e m i n = E min / (dhxL) = —l^jajd and there is a sad- 
dle point in Z at (E)/(dlnL) = e sp — —(3e m i n 2 /2: since 
e sp must be larger than e m i n (as Sl((E)) cannot become 
smaller than 1) the saddle point cannot be valid below 
T c = 1/ f3 c = e rnin /2 = \J a J d and the system freezes in 
low lying states. Although this argument implicitly re- 
lies on using In fl(E) instead of \nfl(E) it does give the 
correct picture for the REM, as shown in ||. 

This picture generalizes to correlated potentials pro- 
vided Ti(r) decreases fast enough at large r. The decay 
must be faster than 1/ In r (which is a rather slow de- 
cay) as indicated by the theorems recalled in App endix 
[b| or also by a simple argument given in Appendix B 2 c. 
Finally, let us point out also that another way to ob- 
tain a transition for SR disorder is to take the d = oo 
limit before taking the large L limit: there the model 
(even without rescaling) always exhibit a transition (in 



G 



the statics and in the dynamics). 



(ii) long range correlations: when the typical V(r) — 
V(r') grows with distance as a power law T(r) ~ \r\ s , 
there is only a low temperature phase and no transition. 
The particle is now always localized near the absolute 
minimum of the potential in the system at r m j n . The 
typical minimum energy Kni n grows as ~ -L s / 2 and thus 
overcomes the entropy ~ Td In L which is never sufficient 
to delocalize the particle. The structure of this single low 
temperature phase is simple: there are no quasi degener- 
ate minima separated by infinite distance (and thus also 
by infinite barriers) in the thermodynamic limit. As can 
be seen on Fig. |3| there is typically a single minimum, 
with many secondary ones near it, but none far away. 
More precisely, as L — > oo, the probability that the low- 
est energy excitation AE(R,L) above the ground state 
(a distance at least R ~ L c from r min ) be smaller than 
a fixed finite (arbitrary) value decays algebraically to 
with L (and AE typ and AE increase algebraically with 
L) . This is the scenario familiar from the droplet picture 
|f, with Proba(AE < T) ~ TL~ S / 2 (i.e in some config- 
urations which become more and more rare as L — > +00, 
there are two far away quasidegenerate ground states). 
In some cases, e.g. in d = 1 Sinai's model (S = 1) the 
distribution of rare events with quasi degenerate min- 
ima has been studied extensively pl|-[43[ . For instance it 
has been shown 0,fl2| that there is a well defined limit 
distribution Q(R)dR[when L — > +00) to find quasi de- 
generate minima E3] at fixed distance between R and 
R + dR, with Q{R)~ R~ 3/2 at large R. 
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FIG. 3. A typical random potential configuration for al- 
gebraically growing correlations 
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FIG. 4. A typical random potential configuration for log- 
arithmic correlations 



Given the growth of correlations one sees that the 
typical energy differences over a distance L scale as 
(V'(O) - V(L)) typ — ±V4ct \nL. Computing the min- 
imum energy is a harder task here, but if one esti- 
mates it as in JlTj through the REM approximation 
l/L d = f_ m ' n P 1 (V)dV (which neglects correlations), 

one finds that it behaves as V m in ~ —2\/ad\nL (for 
Gaussian disorder). This estimate appears rather uncon- 
trolled here since correlations grow with distance, while 
the theorems for uncorrelated random variables apply a 
priori only for correlations decaying slower than 1/lnr. 
In fact the situation is a bit more complex, and as we 
will find below from the RG and our numerics, the lead- 
ing behaviour of V m i n with In L is still correctly given by 
the REM approximation, although the next subleading 
-universal- correction is not. Thus the energy of the min- 
imum — 2v / crdhi L can now balance the entropy of typical 
sites Td In L which yields the possibility of a transition. 
The REM approximation of the model indeed yields a 
transition at T c = \J a jd between a high temperature 
phase for {3 < (3 C = yojd and a frozen phase > j3 c . 
This scenario is confirmed by various approaches in the 
following sections. 

An interesting feature of this model is that the low 
temperature phase exhibits a non trivial structure. Un- 
like long range disorder discussed above, for logarithmic 
correlations we find that the low temperature phase is 
dominated not by one, but by a few states in the thermo- 
dynamic limit. This is in stark contrast with the standard 
droplet picture and is reminiscent of the replica symme- 
try breaking phenomenology, even though we are dealing 
here with a very simple finite dimensional system. 



(iii) marginal case, logarithmic correlations: the most 
interesting case is when correlations grow as T(r) ~ 
4crln|r|. A typical logarithmically correlated landscape 
is illustrated in Fig. |]. One can already see that, con- 
trarily to Fig. ^, it has states with similar energies far 
away. 



One can visualize the transition, and the peculiar na- 
ture of the low temperature phase in Fig. ||,^|, where a 
typical Gibbs measure p(r) is shown in both phases : is 
fairly delocalized at T > T c (Fig. ||) but peaks around a 
few states when T <T C (Fig. ^) separated by a distance 
of the order of the system size. 
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FIG. 5. Gibbs measure in a typical sample in the high tem- 
perature phase (/3 = 0.5 < /3 C = 1), L = 4096 
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FIG. 6. Gibbs measure in a typical sample in the low tem- 
perature phase (ft = 3. > fi c = 1), L = 4096. Only points 
such that pi > 10 -7 are indicated. 

This peculiar nature of the frozen phase can be tested 
by showing that distant secondary local minima with a 
finite AE exist with finite probability in the thermody- 
namic limit. Thus we have investigated numerically the 
distribution Pji^iAE) of the lowest excitation. As il- 
lustrated in (Fig. |J), if the phase is non trivial, we ex- 
pect that this distribution has a well defined limit for e.g. 
R = L/3 when L — * oo with a finite typical AE. Contrar- 
ily to the LR disorder, we expect the probability that e.g. 
AE(L/3, L) be smaller than a fixed number to saturate 
(not to decrease) as L — > oo, i.e that there is a fixed prob- 
ability that a second state within AE exists far away (as 
was already apparent in Fig. We show in Fig. [?], Fig. 
^ and Fig. ^ numerical evidence that this distribution 
has a well defined limit (the details of the simulation are 
discussed in Section fv\ ). Finite size effects are clearly 
important in this system, but their magnitude appears 



compatible with the predictions of our RG approach, as 
discussed below. Thus we conclude that the numerics 
are consistent with the existence of such a limit distri- 
bution and hence with a frozen phase with a non trivial 
structure. 
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FIG. 7. Averaged energy difference AE between the ab- 
solute minimum at r m i n and the constrained secondary min- 
imum (i.e the minimum over the set \r — r m i n \ > L/3), as a 
function of the system size 
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FIG. 9. Distribution of AE in log scale. 



III. RENORMALIZATION GROUP APPROACH 
A. Idea of the method 

We now study the model (|l|, ^ using a renormalization 
approach introduced by us to study d — 2 disordered XY 
models |ll],[lj| . There, one is led to study a neutral collec- 
tion of interacting ±1 charges (XY vortices) in a random 
potential ±V(r) with (^). The single particle problem 
studied here amounts to restrict the Coulomb gas RG of 
|l8| , [l9f to the sector of a single +1 charge. Here however 
there is no charge neutrality and one must be careful to 
study a system of finite size L, as some quantities (such 
as V(r) 2 ) explicitly depend on L, while appropriately de- 
fined quantities have a well defined thermodynamic limit. 

The idea is first to formulate the problem in the con- 
tinuum, with a short distance cutoff a: 

Z = J ^e-W) (16) 

and an appropriately defined cutoff-dependent distribu- 
tion for V{v), and second, by coarse graining infinitesi- 
mally, to relate the problem defined with a cutoff a' = 
ae dl to the problem with a cutoff a. In general, this im- 
plies to be able to follow under this transformation the 
full probability measure of the potential V(r), which is 
quite difficult, as complicated correlations can be gener- 
ated under coarse graining. In some very favorable cases, 
for instance in the d = 1 Sinai landscape (where V(r) per- 
forms a random walk as a function of r - case 8 = 1) , it 
is possible to follow analytically an asymptotically exact 
RG transformation (in the statics and in the dynamics 
p3| ). There a very specific real space decimation pro- 
cedure is required, which can in principle be extended 
here, although it may not be tractable beyond numer- 
ics. The present case of the logarithmically correlated 
potential is thus a priori less favorable but still, thanks 
to some known properties of the Coulomb potential, a 
RG method a la Kosterlitz can be constructed which, we 



argue, should describe correctly all the universal prop- 
erties of the model. There are two possible derivations, 
one which uses replicas and is more precise, and the other 
one without. We start with the latter, which is physically 
more transparent. 

The key observation is that before (and also after) 
coarse graining, the logarithmically correlated disorder 
studied here can naturally be decomposed into two parts 
as: 

V(r) =V>(r)+v(r) (17) 

where V > (r) is a smooth gaussian disorder with the same 
LR correlations as the initial V(r) which represents the 
contribution of the scales larger than the cutoff ae l , and 
v(r) is a local short range random potential which rep- 
resents the contribution of scales smaller than, or of the 
order of, the cutoff ae l . In the starting model v(r) ap- 
pears naturally as a gaussian variable (see below). Af- 
ter coarse graining, v(r) does not remain gaussian, but 
it does remain uncorrelated in space (i.e correlations of 
short range o). The decomposition jl7| ) allows to follow 
the distribution of the V(r) under coarse graining in a 
tractable way. 

The precise way of decomposing the disorder in ( |l7j ) 
depends on the details of the cutoff procedure, but 
should not matter as far as universal properties are con- 
cerned. For illustration let us indicate a simple way to 
do it, a more detailed discussion is given in p9| . It 
starts with the well known continuum approximation 
in d = 2 of the lattice Coulomb potential T(r — r') « 

4d (ln(l^i) + 7 ) (l-^( a )(r-r')) where 6^(r-r') = 

1 for |r — r'| < a and otherwise (7 = ln(2v / 2e c ') and 
C = 0.5772 is the Euler constant). This decomposition 
can be performed more generally, e.g. with other short- 
distance regularization of the potential T(r) (which pre- 
serve the large distance logarithmic behaviour) and in 
any d, which amounts to modify the value of 7. Using 
this approximation the bare disorder (^) can indeed be 
rewritten equivalently as a sum (0) of two gaussian dis- 
order V > (r) and v(r) with no cross correlations and with 
respective correlators: 

(V>(r) - U>(r')) 2 = 4gln * r - ^ (1 - <5 (a) (r - r')) (18) 

a 

v(r)v(r') = 2ajS (a) (r - r') (19) 

With this definition, the problem to be studied is 
rewritten as: 

Z = [*LL z(r) e-^ v> « z(r) = e-' 3v ^ (20) 
J a d 

We can now study the behaviour of the model under 
a change of cut-off. There are two main contributions 
from eliminated short length scales variables. The first 
one can be seen most simply by rewriting the correlator 
in ©: 
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(V>(r) - V>(r')) 2 = 4e7ln ' r * ' (1 - S^'Ur - r')) 

a 

+ Acrdl {\-8 {a '\v-r')) (21) 

explicitly as the sum of a new LR disorder correlator 
with cutoff a' = ae dl and a SR disorder correlator (we 
have discarded terms of order 0(dl 2 )). Thus the origi- 
nal problem with cutoff a can be rewritten as one with 
cutoff a' with (i) a new gaussian LR disorder with iden- 
tical form of the correlator (18) with a replaced by a' 
(ii) a new short range disorder v(r) — > v(r) + dv(r) with 
dv(r)dv(r') = 2(rdlS^(r — r') since it is clear from ( |2l| ) 



that when a 



r,<H 



the LR disorder produces an additive 



gaussian contribution dv to the SR disorder. 

The second contribution resulting from a change of cut- 
off is that neighboring regions will merge. Points rj and 
T2 previously separated as a < |ri — r2 1 < ae dl should now 



be considered as within the same region. The second im- 
portant observation is that the resulting transformation 
can only affect the SR part v(r) of the disorder. Indeed, 
in the region a < |ri — r2 1 < ae dl the LR part V > (r) 
can be considered as constant up to higher order terms 
of order dl. One must view this coarse graining as result- 
ing in a "fusion of local environments" : the two local 
partition sum variables z(ri) and z(r2) combine into a 
single one z(r) according to a rule which we will write as 
z(t) = z(ri) + z(r 2 ). The exact choice of the form of this 
fusion rule is again dependent of the cutoff procedure and 
thus to a large extent arbitrary. 

Putting together these two contributions we obtain the 
following RG equation for the distribution Pi(z) of the 
local disorder z = e~ l3v variable (also called "fugacity" 
in the Coulomb Gas context). 



d t P(z) = p 2 a (1 + zd z f P - dP(z) + d / P(z')P(z")S(z - (z' + z")) 



(22) 



This equation also describes the evolution of the uni- 
versal part of the total free energy distribution with the 
system size. Indeed, the total partition function can be 
written at any scale as: 



Z(f3) 



d d r 



-0V(r) 



d d r 

{ae l ) d 



zi(r)e 



-0V>(r) 



Z( t .) 

(23) 



where the Z;(r) are independent variables distributed 
with Pi(z) and the VJ > (r) are gaussian distributed as 
(|l8|). In the last equality we have coarse grained up to 
the system size : L = a e l .At this scale, there remains a 
single site of (random) fugacity zi* . Thus the distribution 
function of the partition function Z(/3) can be deduced 
from the distribution of the random fugacities at scale I*. 
The distribution of the free energy F = —T\nZ is thus 
given by Pi* (v = F) (where P(v)dv = P(z)dz from the 
change of variable from z to v = —Tin z). Note that the 
~ in ( p3| ) means that these distributions are the same a 
priori only up to subdominant non universal terms (mul- 
tiplicative for Z and additive for InZ). 

For a fixed system size L, the above RG equation de- 
scribes the evolution with the scale I smaller that I* of 
the distribution of z(r), which is the local partition sum 
over scales around r smaller or equal to ae l (i.e of a "lo- 
cal free energy" — Tlnz/(r) = u/(r)). The remaining long 
wavelength disorder at that scale, V^ > (r) should still be 
taken into account when computing the total partition 
sum. 

It is striking that the equation (|3l|) is identical to the 
RG equation for the partition function of a continuum 
version of a directedpolymer on a Cayley tree (a so called 
branching process || ) . We note that it has been derived 
here for a problem with complete (statistical) transla- 
tional invariance, with no ad-hoc assumption about an 



underlying tree structure and simply adapting to the 
present problem the Coulomb gas renormalization a la 
Kosterlitz. That the correspondence between the two 
problems naturally appears within the RG with no addi- 
tional assumptions, is even more apparent on the deriva- 
tion using replicas of the next section. Thus we consider 
that this establishes on a firm footing the strong connec- 
tion between the two problems. 

Before analyzing the consequences of the above RG 
equation let us sketch the more precise derivation using 
replicas. Other derivations without replicas are also pos- 
sible and we refer the reader to |l^] for more details. 



B. derivation of the RG equation using replicas 

Let us consider the whole set of moments Z m which 
encode for the distribution function P\Z\. They can be 
written as: 



d d v x d d r m 



a" a" 
This can be rewritten as: 



f d d Y\ d d r„ 

Z m = 



(E <=1> .. m v (ri )) 2 (24 ) 



e 4 



£«* J =i...«r(r ( -r J ) e mW'ln£ (35) 



We have used that T(r, r) = T L (0) = 2<rln^. One can 
choose a regularisation, e.g. T(r — r') = V(r)V(r') — 

— crln 1 ~ J} . Notice that only the large distance be- 
haviour of the above correlator is important for the fol- 
lowing renormalization. 



10 



We now switch to another representation of the replica 
partition sum. (|2^) is a partition sum of m particles 
located at ri,..r m corresponding to m replicas. Now 
instead we will index the configurations using (vector) 
columnar replicated charges. To each point r, within 
a hard core size a, we associate a m-component vector 
n whose components n l (r) are either 1 or depending 
on whether the particle corresponding to the i-th replica 
is present within a of r (|r — r^| < a) or not. These 
charges thus correspond to n = (0, 1, 0, ..0, 1, 1) since sev- 
eral replicas can be present near a given point. Choosing 
a columnar hard core for the vector charges corresponds 
to a choice of cutoff, which is arbitrary, but the universal 
features of the renormalization should not depend on it 

H- 

The m— th moment of P[Z] then read 



L 



2 am 2 I 



{™Ll Q 



(26) 



d d r n 



|r Q -r Q ,|>a a 



■ exp 



-2/3 2 cr n a n a i In ( — 



a<a' 



where the primed sum correspond to a sum over all dis- 
tinct non zero configurations of replica charges n a at sites 
r Q . We have defined n a — ^ n z a as the total number of 
replicas present in a given charge (n' l a = 1). The quan- 
tities y[n] are functions of the local vector charge and 
are the so-called vector charge fugacities. In the bare 
model they appear as soon as the continuum approxi- 
mation to the lattice Green function is used and read 
Y[n] = e -2<T7 " . Since we are studying a single parti- 
cle problem, there is also an important global constraint 
on the configuration sum that only one particle in any 
replica i is present in the system, i.e: 



i 



(27) 



which is preserved by the RG. 

The RG equations for this model read: 



«9,y[n] = (d + f3 2 an 2 )Y[n] 



(28) 



Sd- 



i J2 Y[n']Y[n"] 



n'+n"— n 



where the sum is over n' and n" non zero vector charges 
(also n is non zero) and Sd-i is the volume of the unit 
sphere in dimension d. We recall that n = ^i=i ri i- 
These equations are obtained by a generalization of the 
Kosterlitz procedure |46) as follows. The first term comes 
from explicit cutoff dependence in ( |26] ) . Upon increasing 
the cutoff infinitesimally a — > a' = ae dl the integration 
measure and the a dependence in all logarithms combine 
to give Y[n a ] -> Y\nje d ^ d+a02 <y. We have used that 



2 J2 a <a' n a"n a ' = m 2 - J2 a n a which holds due to (p7[). 
The last term in the above equation (08) comes from the 



fusion of replica charges upon increase of the cutoff. The 
above RG equations hold for any m. 

We should now look for solutions of this set of equa- 
tions analytically continued to m — * 0. One way to do 
that is to find a convenient parametrization for the set 
of Y[n]. Here we preserve replica permutation symmetry 
within the RG and we can thus choose Y[n] to be a func- 
tion of n = rii only. Then we define the parametriza- 
tion Y[n] = J dz$i(z)z n = J du^i(v)e~ Pnv . The differ- 
ent terms in the above equation then translate into 



n 2 Y[n] = J dve-^^dvf&iiv) (29) 



Y, Y[m]Y[n 2 ] 



z") 



-2Af$i(z) + S(z)JV 2 (30) 

where M = J One then easily converts the 

equation for <$>i(z) into an equation for a normed func- 
tion Pi(z) = $>i(z)/Af> defined only for z > 0, with 
W> = f z>0 ®i{z) (as in Q) by noting that W> con- 
verges quickly to A/> = 2d/Sd-i- The resulting equation 
for Pi (z) is exactly the one (^2|) given above, and its phys- 
ical interpretation in terms of the probability distribution 
of the fugacity (i.e the local partition sum) was given in 
the previous section. 

What is the small parameter which controls the valid- 
ity of the above RG equations (with and without repli- 
cas)? In conventional Coulomb gas context, these RG 
equations are known to become exact in the dilute limit 
of non zero (vector) charges |l6| . It is easy to see that this 
corresponds to the tail of the distribution P(z) for large 
z (or equivalently small v). This is further confirmed, 
a posteriori, by the remarkable universality properties 
of the resulting non linear RG equation (|22|), analyzed in 
the following section, which arises precisely in this region 
of z. So to obtain the universal behaviour (e.g. of the 
distribution of free energy) we are working with sufficient 
accuracy. On the other hand the bulk of the distribution 
Pi(z) seem to be sensitive to details of the cutoff pro- 
cedure (e.g. details in the fusion rule) and as discussed 
below, is thus likely (unless proven otherwise) to be non 
universal. 



C. analysis of RG equation and results 



1. KPP front propagation equation and velocity selection 



Let us analyze the solutions to the RG equation (|22 
In terms of the (local free) energy variable v(r) — 
— Tlnz(r) (from ( p0| ) and its distribution Pi(v) = Pi(z 



-0v 



)0e 



-0V\ 



it has a well defined zero temperature limit, 



since then the fusion rule simply becomes the extremal 
rule v' = min(w 1 ,w 2 ) leading to : 



d z P(v) = ad 2 v P + dP(v) -1 + 2 



+00 



P(y')dv' (31) 
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To be able to work at all temperatures, it is in fact useful 
to trade the distributions Pi(z) or Pi(v) for the generat- 
ing function @j47j]: 



Gi.p(x) = (e ze ) Pl{z) = (e" 



>Pl(v) 



(32) 



We will sometimes drop the index (3. At zero temper- 
ature, the double exponential becomes a theta function 
and Gi(x) simply identifies with the distribution func- 
tion: 



Gi-p =+ca (x) 



+ OC 



Pi(v)dv = Proba(u > x) (33) 



and for all /3 it is a decreasing function of x with 
Gi(x — > — oo) = 1 and Gdx — ► +oo) = 0. Note 
the asymptotic behaviour |4§f| at very large negative x, 
1 — Gi(x) ^< z >p, e^ x . The temperature appears only 
via the initial condition || and the problem at hand is 
thus to determine the large I behaviour of Gi(x) for a 
given initial condition. 

The equation ( |3l| ) is easily transformed, at all temper- 
atures, into the Kolmogorov (KPP) non linear equation 



^diG{x) = ^d 2 x G + F[G] 
F[G] = -G(l - G) 



(34) 
(35) 



which describes the diffusive invasion of a stable state 
G = into an unstable one G — 1. This class of 
equations admits a family of traveling wave solutions 
Gi(x) = g(x + m(l)) which describe a front moving to- 
wards negative x and located around x ~ —m(l). This is 
readily seen by plugging this form in ( |34| ) and assuming 
that dimp(l) — > c one obtains the equation for the front 
shape: 



l -cg'{x) = ^g"(x)+F[g(x)] 



(36) 



The family of such traveling wave solutions g c {x) can 
thus be parametrized by the velocity c. (|36| ) simplifies 
for large negative x when g ~ 1. Denoting g = 1 — g and 
using that F[g] ~ — g for g w 1, one finds the linearized 
front equation for g c : 



1 cr_„ 
d Cff = d 9 
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(37) 



This equation allows to relate the speed of the front c to 



the asymptotic decay of the front, since if g{x) 
large negative x one finds: 



1 

Q 



for 



(38) 



The problem at hand now is to determine toward which 
of these front solutions g c (x) will Gi{x) converge at large 
I, and thus what will be the asymptotic front velocity. 
This velocity will determine the intensive free energy of 
the original problem. Indeed, the convergence at large I 



of the solutions of non linear equations of the type ( |34| ) 
(with a general -F[G]) towards one of such front solu- 
tions, and the corresponding problem of the selection of 
the front velocity c, is a famous problem, still under cur- 
rent interest in nonlinear physics ^-|53|. 

The simplest argument is to use the fact that for very 
large negative x, one must have g(x) ~ e@ x and thus 
a = /3. This seems to imply that the front velocity is: 



(39) 



This however is not always true. First note that the 
curve c(/3) has two branches, i.e that in this naive es- 
timate two different (3 would correspond to the same 
velocity. The special point C = \fd~fa corresponds to 
c = c* = 2dy/ a/d. For more general non linear equa- 
tions one usually relies on the so called marginal stabil- 
ity criterion (e.g. which shows that the large f3 branch 
is unstable and can be eliminated) Here there are 

rigorous results available : the Bramson theorem ]54| en- 
sures the following results, which are independent of the 
precise form of F[G] (up to some rather weak conditions 
on F[G] ||): 

(i) At high temperature, (3 < (3 C = \J dja the asymp- 
totic front is indeed an exponential for large negative 
x and Gi(x) uniformly converges towards the traveling 
wave solution g c (i3){x + m(l)) where the velocity is given 
by ([59|), thus continuously dependent on temperature. 

(ii) At low temperature (3 > [3 C the velocity freezes to 
the value c — c* and the front decays as: 



-xe 



(40) 



for large negative x, thus independent of the tempera- 
ture. The solution G\{x) uniformly converges towards 
the traveling wave solution g c *(x + m(l)). Thus in that 
regime, one must then distinguish two regions in Gi{x) 
at large I, the front region and the region very far ahead 
of the front (x + m(l) ^> \fl) where the decay is again 
as Gi(x) ~ exp(/3x) as it should: this will be discussed 
again below. 

There are additional rigorous results from | p4j and in 
particular the remarkable fact that not only the veloc- 
ity but also the corrections to the velocity are universal 
(independent of F[G}) i.e one has for the position of the 
traveling wave mp(l) at "time" I: 



i(l)= (^p + fi-^dl + Cst (3<(3 C 

f3>(3 c 




- { 2dl --Inl 
d\ 2 



- ( 2dl --Inl 
d\ 2 



(41a) 
(41b) 
(41c) 
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2. Results for the fugacity and free energy distribution and 
extremal statistics 



These results on the KPP equation ( |34| ) can now be 
translated (via (|3^)) into results for the fugacity distri- 
bution Pi(z) and for the distribution of free energy Pi{v). 
One finds that Pi(z) and Pi(v) also take the form of a 
front at large I, e.g.: 

Pi(v) ^p(v + m(l)) (42) 

with p(v') related to g{x) by g{x) — J ,p(v')e~ 
Thus we obtain that the local free energy is: 

- /T 1 (In z) ~ -mp{l) (43) 



can be confident that the RG approach gives the exact 
result) and what may not be (as it depends on the details 
of the cutoff procedure, yielding e.g. a different KPP non 
linearity F [G]). 

Let us start with T = 0. We find (cf. (flj,[45c|)) that the 
minimum V m in of L d logarithmically correlated variables 
behaves as: 

Vmin = -2VadhiL + ^J^InQnL) + SV (46) 

and SV is a fluctuating part of order 0(1) ■ Since at T = 
one has p(v) = g'(v), from the result ( |40| ) we get that the 
tail of the distribution of u = SV— < SV > for u — ► — oo 
is universal and behaves as: 



up to a finite constant, w here the position of the front 
iri)3(l) is given above in (41a). Using the result (23), 
TV = d\n(L/a) = dl* we obtain using the RG that the 
free energy F[V] of the system of size L reads: 



p(u) 



(47) 



F[V] =h W In L + SF 



(44) 



where SF is a fluctuating part of O(l) of probability dis- 
tribution p(SF) and the intensive free energy reads: 



h{P) = 
h(P) = 
h(P) = 



+0 



1 

hTT 



/3 C V 2 dlnL 
_1_ f 3 1n(lnL) 
/3 C V 2 dlnL 



O 



O 



P<Pc = \l- 



SI' ^ 



(45a) 
(45b) 
(45c) 



with f3 c = \J dj a . Thus we find a distribution different 
from the Gumbell distribution, and thus correlations do 
matter. 

The question of what is universal in this distribution 
is non trivial. We find from our method that the full 
distribution of P(u) depends on the detailed form of the 
front (and thus on F[G] and a priori on the cutoff proce- 
dure) and is thus less likely to be universal (although this 
remains to be investigated). Hence we believe that uni- 
versal features include at least the tail of the distribution 
©• 

The above result ( |47| ) carries through the tail of the dis- 
tribution of the free energy u = F— < F > for u — > — oo 
for T < T c and it was shown in || that for T > T c one 
has: 



where the factors 1/2 and 3/2 which arise in the finite 
size corrections are universal. 

Thus we have found using our RG method that in any 
dimension d > 1 the original model (|l]||) exhibits a phase 
transition at (3 = f3 c (d). This transition is very similar 
to the freezing transition of the continuous version of the 
random directed polymer on the Cayley tree. Our RG 
thus confirms that the REM approximation ( |lO| ) to the 
model does give the transition at the same (3 Cl and with 



p{u) 



(3<(3 C 



(48) 



same asymptotic intensive free energies (lib) as (45c). 



It allows however for a more detailed study and shows 
that the universal finite size corrections differ in the two 
model. In the REM the above formula with the factor 
1/2 holds in all the low temperature phase, which is not 
the case for the present model. Thus the present model 
is in a different universality class than the REM. The 
physics that we find here is much closer to the one of the 
directed polymer on the Cayley tree: it remains to be 
seen whether this can be extended to other observables. 

The RG method also yields the distribution of the O(l) 
fluctuating part SF of the free energy, and in particular 
at T = it gives a result for the extremal statistics of the 
correlated variables. We must now carefully distinguish 
between what is clearly universal (and thus for which we 



D. More on fronts, REM via nonlinear RG and 
extremal statistics 

To illustrate how the previous results fit in a broader 
context, let us show how the simpler properties of ex- 
tremal statistics of uncorrelated variables and of the ran- 
dom energy model can be recovered within the same RG 
framework. This provides, en passant, yet another solu- 
tion of the REM. 



1. uncorrelated variables with fixed distribution: Gumbell via 

RG 

Let us consider N = e ld = (L/a) d independent random 
variables V(r) r = 1, ..N with a fixed distribution P(V) 
(d here does not play any role as the true variable is Id 
but we keep it for the sake of comparison) . The generat- 
ing function of the distribution of the partition function 
Z \V} = Er e-^M of model (@) reads: 
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G t {x) =< eM~Z[V]e^) > P{V) 
J dVP(V)exp(-e^ x - v ^ 



It satisfies the equation: 



^lnlni = l 



(49) 



(50) 



Or, interestingly enough, it obeys a KPP type equation 
with no diffusion term: 



1 



diG = F[G] 



F[G] = G\nG 



(51) 
(52) 



The Gumbell distribution now emerges naturally from 
the front solutions of this equation. Writing Gi(x) ~ 
g(ai{x + mi)) and assuming di(aimi) — > c yields eg 1 = 
<7 In g whose solutions with the above boundary condi- 
tions are g(y) = exp(— r ye v ^ c ) (7 being a positive con- 
stant). We have assumed diai 0. Since there is 
some freedom of choice for ai and mi, one can always 
set c = 7 = 1. The determination of the rescaling factors 
ai and m(l) is performed in Appendix |^. At T — one 
has P(V m i n ) = —G'(V m i n ) and one recovers the known 
results from probability theory for the convergence to 
the Gumbell distribution detailed in Appendix |^, but 
the generating function Gi (x) takes a Gumbell form also 
at finite T. 



2. REM via RG 

We now turn to an alternative derivation of the solu- 
tion to the Gaussian REM model using a RG approach 
and a traveling wave analysis. This allows to make some 
connections with the correlated case studied previously. 
Let I — InL and lnA^ = Id. 

We want to write a RG equation for: 



G,(x) 



Pl(V) 



(53) 



where the single site distribution Pi (V) now is scaled with 
I. We introduce 

Gi(x) = (^ X ~ V) ) PiV =exp(e- w lnG i (x)) (54) 

Let us choose the single site distribution Pi(V) which cor- 
responds to the REM approximation ( |To| ) of the model 
studied here i.e the gaussian: 



Pi{V) 



1 



-e 4<7i 



It satisfies: 



8iPi(V) = ad 2 v Pi{V) 



(55) 



(56) 



(57) 



One easily checks that it implies that: 

diGi{x) =ad 2 x Gi(x) 
This leads to the equation for Gi(x): 

diG = ad^G + dG\nG - a(l - e- ld )^{d x G) 2 (58) 

Thus the RG equation of the REM, for large I reads: 

diG = ad 2 x G + dG\nG - a^{d x Gf (59) 

and is almost a KPP equation, except that it has an ad- 
ditional gradient (KPZ type) term. This term here plays 
an important role and yields a different universality class 
than KPP. We now search for the front solutions. 

Let us rewrite the exact equation (p8|) using the func- 
tion h = — InG (remember that < G < 1): 



drh = dh + ati' 



ae- ld h' 2 



(60) 



For large I we can neglect the decaying nonlinear part, 
and we now look for a solution of the linear equation. 
The only front solution of the form h(x) = h(x + rn(l)) 
with dim(l) — > c which satisfies the boundary conditions 
h(— 00) = and h(+oo) = +00 is the exponential: 



hi(x) = e o(l+m(!)) 

dim(l) = c = hi 

a 

fix 



(61) 
(62) 



By using again the hi(x) ~ e^ x boundary condition at 



-00, we find a — f3 and: 



c(/9) = g + <r/3 



(63) 



as in (p9[). This is correct in the high T phase and 
yields the correct REM value for the in tensi ve free en- 
ergy /(/?) = c(/3)/d + 0(1/ InL) as in ( JTTb| ) (and also 
correctly yields the absence of non trivial finite size cor- 
rections). Thus for the REM in the high T phase we 
find: 



Gi(x) w exp(-e l3( - x+m{l) '>) 



(64) 



thus again a Gumbell form, with ai = (3 and m(l) = 

To see the transition to a low T phase for (3 > (3 C = 
\Jdja and the freezing of the velocity at c = c* = 2^/da, 
one needs to carry a slightly more detailed analysis (dis- 
carding again the decaying nonlinear part). The general 
solution of the linear part of the equation (pOh is: 



hi(x) 



dx'- 



1 



Jd- 



ho(x') 



(65) 
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where ho(x') can be interpreted as the hi(x') at earlier 
time Iq such that the nonlinear terms can already be ne- 
glected and decays as ho(x') ~ e@ x for x' — > — oo. 

This formula nicely exhibits the REM transition. In 
the high T phase, using the asymptotic form ho(x') ~ 
e /3a; we g nc L that there is a saddle point at x' = x + 2a 01 



This gives h t (x) ~ e P(.*+cW) w ith c(0) given in Q. 
The front is centered at x* = —c(j3)l and consis- 

tency requires that the corresponding saddle point x 1 * 
moves to — oo so that the asymptotic form of Hq(x') can 
indeed be used. Hence we havex'* ~ (a0-£)l. Thus the 
saddle point become inconsistent and the high T solution 
ceases to hold, for > C = \fdj~a. 

The solution in the low T phase is easy to find. Setting 
x = —m(l) + y one finds for large I: 



hi(y) 



Jd- 



■^O'-iMMeiv / dx'e-£ x 'h (x') (66) 



where we have denoted c* = lim;^ +oc m(l)/l and ne- 
glected the additional factor e~ x /( 4cri ) in the integral. 
This is correct provided the integral: 



dx'e " 2 " x ho(x') 



(67) 



is convergent, i.e c* < 20a. The consistent choice for c* 
and m(l) must be: 



3. conclusion on RG fronts and extremal statistics 

Thus we have seen on two examples that extremal 
statistics problems (and their T > thermodynamic 
model counterpart) can be studied using non linear RG 
equation with traveling wave solutions. In one example 
(uncorrelated rescaled variables, i.e the REM) the RG 
equation is exact, while in the second (logarithmically 
correlated variables) we only know it presumably in the 
tails. The front position represents the typical value of 
the minimum Vmin as a function of I = d~ 1 hiN while the 
shape of the front gives the distribution of the V m in (resp. 
of the free energy F) . This suggests that a broader class 
of such models can be approached by these methods, and 
raises the question of universality. 

Studies of such non linear equations [|53| usually dis- 
tinguish between pushed fronts where the velocity relaxes 
exponentially in I (velocity selection by non linear terms) 
and pulled fronts (velocity selection by the marginal sta- 
bility criterion). The extremal statistics (and the glassy 
phase) correspond to the pulled fronts. There one ex- 
pects a very broad universality as stressed in |5l],|52j : not 
only is the asymptotic front universal but also the ve- 
locity and its corrections. In a nutshell, the argument 
for the universal | In I corrections to the front position 
comes from matching of the universal tail of the front 
g(y) ~ [Ay + B)e~ l3aS - with y = x + m(l) with the far tail 
region, so far ahead of the front that one can linearize 
the KPP equation and get: 



2Vad 



m(l) = Uld - ~ ln(47rcr/) j + (H I ) ((is i 



1 - Gi(x) « e-^V(y) diMv) = °dlMv) (72) 



which ensures that ( p6|) has a proper limit hi(y) ~ 

Ae^y = Ae^ v which is again a Gumbell form for Gi(x) 
but now is temperature independent. This holds for 

V > p c = VW^- 

From this method of solving the REM we have recov- 
ered the result of || namely that for (3 > (3 C the free 
energy behaves as: 



fL{P > ~ dl ~ fa V 2 dhxL 



The only matching solution is ipi(y) = y/l 



(69) 



In addition we recover, for T = the result for the min- 
imum V m in in the REM approximation: 



V min = -2VadlnL+ -W^ln(lni) 



SV 



(70) 



with u = SV— < 5V > distributed with a Gumbell dis- 
tribution: 



Proba(u > x) = exp(-Ae /3cX ) 
where A is a constant. 



(71) 



_ ; ; 3/2 e -y 2 /(4al) _ 

Inserting y = x + mil) = x + c*l + ClnZ immediately 
yields C = 3/2 for proper matching. As discussed in ]5^j 
this universality extend for pulled fronts in a very broad 
class of non linear (or coupled non linear) equations and 
holds for steep enough initial conditions (i.e in the glass 
phase in our language). 

This argument fails in some cases, such as at the bifur- 
cation between pushed and pulled fronts (e.g. at the glass 
transition = C or equivalently when the initial condi- 
tion has slow decay ~ exp(— c x)) (see e.g analysis in 
pof). Interestingly, it clearly fails also for the non linear 
equation corresponding to the REM model, which is thus 
in a different universality class (this may be related to the 
fact that fronts are here unbounded |5q1 ). Presumably 
what happens there is that the coefficient A vanishes, and 
the solution is exactly e~^ ay , hence the ^lnZ (since the 

above matching function is now tp(y) ~ l~ 1 / 2 e~ y2 ^/( 4o "')). 

Next is the question of universality. We will address it 
only for our model of gaussian variables with logarithmic 
correlations. We have recast the RG equation (g2[) into a 
KPP equation with a specific non linear term F[G] . From 
our RG we have obtained F[G] = —G(l — G). The struc- 
ture of the RG derivation suggests that we have obtained 
correctly the two lowest orders of F[G}. From the above 
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discussion this is enough for the universality. Thus, and 
we call it the restricted universality scenario, it is likely 
that higher order terms F[G] = -G(l -G)+0((1- G) 3 ) 
are non universal and thus that only the tail of the dis- 
tribution of the minimum of log-correlated variables is 
universal. 

Let us mention however that we were not able to rule 
out another scenario, the broad universality scenario, 
such that the true distribution of the minimum of log- 
correlated variables is indeed universal. If this was true 
the following conjecture would be tempting: since we 
know that for uncorrelated variables the KPP RG equa- 
tion is exact with F[G] — GlnG and a = (and see 
Appendix [b] is asymptotically exact even for weakly cor- 
related ones), one could conjecture an interpolating KPP 
equation ( p4| ) with F[G] — GlnG and a > which would 
gives exactly the distribution of the minimum of log- 
correlated variables. Unfortunately we have been unable 
to confirm (but also to strictly rule out) numerically this 
conjecture, due to the very large finite size corrections, 
as discussed in Section IV. 



The RG analysis performed in this paper indicates that 
adding correlations will not change this conclusion. In- 
deed, one first coarse grains up to scale lo = hi(L) — g In 2. 
At this scale, the system can be described by two local en- 
ergies (one for each half) of minima v\ and u 2 distributed 
according to Pi (v), to which should be added a term SV 
which correlates the two halfes and is gaussian of vari- 
ance ~ ^ln2. This however does not change the fact 
that the difference V m in\ — V m i n 2 ~ 0(1). Thus one still 
finds that there exist secondary minima of O(l) in energy 
from the minimum, and a typical distance L a way f rom 
the absolute minimum. As discussed in Section II D this 



property was also confirmed by numerical simulations. 

It is natural, in view of the analogy with the directed 
polymer on the Cayley tree, to introduce the "overlap" 
between two different states (i.e positions of particles) ri 
and r 2 as: 



<j(ri,r 2 ) = 1 



ln(a + |n - r 2 |) 
hxL 



(73) 



We expect it to be non self averaging and characterized 
by the "overlap distribution" : 



E. structure of low temperature phase and replica 
symmetry breaking 



(74) 



Let us now return to the structure of the low tem- 
perature phase for the particle in the d-dimensional ran- 
dom potential with logarithmic correlations. We argue 
that (i) it has a non trivial structure, with a few states 
(ii) this structure is reminiscent of the so-called "replica 
symmetry breaking" |E| This non trivial structure can 
be caracterised more precisely here as the various states 
of the model correspond to the different positions of the 
particle, and have thus a natural meaning in real space. 
In particular, the minima of the "energy landscape" (or 
metastable states) are nothing but the local minima in 
the sample of the random potential for our problem. A 
precise caracterisation of these 'local minima' is given be- 
low. Also, approximate replica solutions of our model are 
shown in the following to exhibit RSB at low T. 



1. Spatial distribution of secondary minima 

Let us start with a simple argument: for a given re- 
alization of disorder, we divide our system into two sub- 
systems of size L d /2, and call V m i n \ and Vmini the two 
corresponding minima in each subsystem. 

Within the REM approximation, we know from 
that V min i-V min 2 ~ (y\-V2)\foJd ~ O(l) where yi and 
y 2 have independent Gumbell distributions. Thus clearly 
in that case there is a non trivial structure: the secondary 
minimum (defined as being constrained to lie within the 
other subsystem) is typically within AE — 0(1) in en- 
ergy of the absolute minimum (and within this approxi- 
mation the distribution is also easily computed). 



Although we have not attempted to compute this func- 
tion directly using our RG it is natural to expect that, as 
in the REM and the DPCT, it is non trivial for T <T C 
and reads: 



P2(Q) = Y <%) + (!- J )S(l-q) 



(75) 



Similarly one expects that in a given disorder environ- 
ment, the probability of finding an overlap q between 
two thermal realizations becomes in the large L limit: 



Y(q)dq=(l-Y)6(q)+YS(q-l) 



(76) 



with Y = 1 — — and Y has the same distribution as in 
the REM. Thus the natural expectation, from the DPCT 
analogy, is that the overlap in the low T phase will be 
either 1 or (i.e secondary minima - of energy differ- 
ence of order T - will be either near the absolute one 
lnri 2 /lnL — > 0, or a distance n 2 ~ 0{L) typically a 
fraction of the system size away) It would however be 
of interest to investigate further these properties in the 
present model, in particular to obtain more detailed in- 
formation at intermediate scales, e.g. correlations prob- 
ing the whole range lnria ~ (ln_L) a with < a < 1. 



2. Approximate replica symetry breaking solutions of the 
model 

Let us now turn to the replica representation and dis- 
cuss how the present model exhibits a form of "replica 
symmetry breaking" . The replicated partition sum reads: 
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d d r! d d r.„ 



-a^'E.^in^V^e™ 2 ^ 21 "^ (77) 



It turns out that various approximations of this partition 
function (specifically the REM and the DPCT approxi- 
mations) are dominated, in the limit m — ► by replica 
symmetry breaking configurations. 

In the context of 2d Dirac fermions with random vector 
potential (see Section VI) an estimate of ([77]) was given 
in [p4| . For small 3 it is clear that the exponential con- 
taining the logarithmic attraction between replicas does 
not decay fast enough and thus the integral is dominated 
by the configurations where the replicas are all far away 
O(L) apart thus: 



md(l+f/3 2 ) 



(78) 



This estimate of Ref. [|24| is in fact incorrect as it misses 
the glass transition. Indeed, one can redo this argument 
using configurations where m/p packets of p replicas are 
0(L) far apart (while in each packet the replica (inde- 
pendent particles) are close to each other). This estimate 
was performed in Ref. |57fl and gives instead: 



Z" 



j3 am' -+-d——{3 2 crrn('m—p) 



(79) 



The interaction term being proportional to the num- 
ber of pairs of replicas in different packets, which is 
m(m — p)/2. In the limit m — > one can then optimise 
over < u = p < 1, i.e: 



Z" - exp I din- max ( - + -B 2 u 

a o<m<i \u a 



(80) 



For 8 < 8 C = y ~ the saddle is for p = 1 and one recov- 
ers the above expression. For 8 > 8 C — — one finds that 
the saddle is for u = B c /8 — T/T c which gives: 



d\n ±2 



(81) 



Thus this calculation yields a transition. In Ref. |57j it 
was claimed that it does not correspond to replica sym- 
metry breaking. We believe that this is incorrect and 
that this is a (one step) RSB estimate of the above par- 
tition sum. This is clear since this calculation exactly 
amounts to the corresponding one for the REM approxi- 
mation of the model, i.e replacing in ( |78| ) Yli<j m J 3 ^ 
by X)i<j (l — ^n^j) ln(L/a). In the REM we know from 
Ref. || that the correct solution for T < T c can be 
obtained by performing the analytical continuation to 
Tji — ► on a RSB saddle point (note that the REM finite 
size correction i In In L is also obtained from the saddle 
integration) . 

One can go one step further and use an argument 
based on universality, which puts the present problem in 



the DPCT universality class (for some observables such 
as the free energy distribution). For the DPCT, it was 
shown in Ref. J10| that one can also recover the correct 
result for the averaged free energy by considering directed 
polymer configurations which break replica symmetry as 
77i — ► 0. It remains to be demonstrated how to obtain 
other universal quantities, e.g. the |lnlnL finite size 
corrections, via a RSB saddle point calculation. 

It is interesting to see how the features associated to 
RSB arise from the RG developed here, despite the fact 
that it is explicitly replica symmetric. Quite generally, 
if one can find independent local free energy variables 
with an exponential distribution P(f) ~ one natu- 
rally obtains a RSB picture. This is the case here, up to 
some more detailed universal preexponential structure in 
P (/). The important feature of our RG is thus that it 
follows the full distribution Pi(z) of local disorder (i.e. of 
local Boltzman weights z) which becomes algebraically 
broad as I — ► +oo. Here this property is sufficient to 
show that the low T phase has a structure reminiscent of 
RSB. Indeed, let us again coarse grain the system up to 
an already large scale L$ = ae l ° but still much smaller 
than L, the ratio L/Lq — e' 1 = M being large but fixed as 
L — > +oo. Assuming that L is so large that Pi (z) has 
reached its fixed point already (except in a remote tail 
region corresponding to very rare events). Since one has 
the decomposition (0^) the RG tells us that the sample is 
divided in M subsystems with free energies Fi — Vi + , 
i = l,..M where the variables z» = x,..m = e ~^ Vi are in- 
dependently drawn from the common distribution Pi (z) 
and the V{" are still correlated but gaussian. Neglecting 
first the V^ > we are left with a system of M subsystems 
of Gibbs measure: 

(82) 



E 



Since the Zi are drawn from a distribution with algebraic 
tails P(z) - l/z 1+ ^ with fi = T/T c one has < z >= +oo 
for T < T c and, as is well known, the partition sum 
is dominated by a few of the zi variables P,|56[| (which 
in essence is the physics associated to RSB). Since the 
correlated V > variables are in finite number and with 
gaussian tails they cannot change the exponential tails 
of the Fi and thus adding them back should not change 
the above conclusions. 

Thus here, although the RG is replica symmetric, since 
it allows for generation of broad tails it can capture fea- 
tures usually associated with RSB. 



IV. NUMERICAL STUDY 

Since we found via the RG and other arguments that 
there should be a transition in any dimension d > 1 it is 
particularly convenient to perform numerical simulations 
in the "extreme case" of d = 1 (i.e the further away from 
mean field). However, even in d = 1 numerical simula- 
tions are delicate because the finite size corrections are 
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very large (and interesting to study, in order to distin- 
guish various universality class). Indeed we have found 
that the main numerical uncertainties come from the fi- 
nite size effects and not from the number of averages. In 
most of the numerical work averaging over ~ 10 real- 
izations of disorder was sufficient, while a simulation of a 
system of size 2 21 ~ 2 10 6 leads to important corrections 
to the thermodynamic behaviour of the model. In view of 
this, we believe that the previous numerical investigation 
p| was at best approximate. 

We have considered a lattice model in d = 1 with 
L = 2™ sites. The potential V(r) on each site (r — 
1, ..L) was computed from its Fourier components V(r) = 

wl/2(~ 1)'' + zJfc=i Wk cos(2nkr / L — (f>k), eliminating 
the k = mode, with Wk independent gaussian variables 
WfcWfc' = A(k)6k,k' (k,k' = l,..L/2) and each <fik inde- 
pendently distributed uniformly in [0, 2ir}. We choose 
A(fc) such that: 



T(r- 



V(r)V(r') 
27Ty? cos(^(r - r')) 
L fe =i |sin(^)|J6-2cos( 



fc=i |sin(f ) |^6-2 cos(^) 



(83) 
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FIG. 10. Zero temperature limit: finite size corrections to 
the minimal energy. Plotted is 2l—\V m in | versus In Z (Z = In N) 



so that (V(r) - V{r')) 2 = 4crln(r - r') forl«r-f / C 
i/2. This is the choice which also corresponds to corre- 
lations along the axis y — on a 2d square lattice. 

The behaviour of the model has been studied, without 
loss of generality, at zero and at finite temperature for a 
disorder strength a = 1 (others value of a can be incorpo- 
rated in the definition of the temperature scale). We have 
first computed the average minimum e m i n = V m i n /hiN 
(with N = L) for system sizes ranging from L = 2 7 = 128 
to L = 2 21 ~ 10 6 and for each size we have taken the 
average over 10 4 realizations of disorder. An estimate 
of the uncertainty on the disorder average was made by 
measuring the variance of a series of average over 10 re- 
alizations. This variance was found to be of the order of 
10~ 3 for all the value of V m i n . The results are plotted 
in Fig. [l0| We recall that the RG prediction reads for 
(7 = 1: 



Next, we have plotted the distribution of V m in in Fig- 
|TT| and compared with the prediction of the RG for the 
tails. Here also the agreement is satisfactory. 
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(84) 



We should first note that if one does not assume any- 
thing about the finite size corrections, the resulting uncer- 
tainty on the ratio e m i n = V m i n /lnN is very large even 
for sizes L — 2 since the ratio § ln(ln N) / In N w 0.3. 
Hence with no assumption it is hard to estimate e m ;„ to 
better than 10 per cent accuracy. 

However, if one assumes that e m % n = 2, the plot in 
Fig. |Io| shows the existence of the ln(lniV) corrections 
with a slope definitely larger than 1 and consistent with 
3/2 (although the accuracy is not excellent). It is how- 
ever sufficient to rule out a REM type behaviour and is 
consistent with the RG prediction 
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FIG. 11. Distribution of E„ 



Fin ally, w e have plotted the "glass order parameter" 
Y2 = "YlirVr which is non zero when the system is dom- 
inated by a few states. It is consistent with a very slow 
convergence towards Y% = (l — T/T c )6(T c — T) but clearly 
other forms cannot be ruled out. 
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FIG. 12. Plot of Yi — X], Pr as a function of temperature 
for different system sizes L — 2 N 



V. RELATIONS WITH LIOUVILLE AND 
SINHGORDON MODELS 

In this Section we describe the relation between the 
problem of the particle in the log-correlated random po- 
tential and the Liouville and sinh-Gordon models. Exact 
results on the sinh-gordon model are compatible with 
(and also point out towards) the existence of the transi- 
tion at (3 = f3 c . 



A. Relations with the sinh-Gordon model in d = 2 

and d — 1 

Let us start with the correspondence with the sinh- 
Gordon model. Although less direct, it is also simpler 
to analyze, as the model does not contain subtle bound- 
ary conditions problems. The interest of the connexion 
is that the sinh-Gordon model is intcgrable in d = 2 and 
d = 1 (Boundary sinh-Gordon) [j58[-|60| . 

The connexion requires introducing a slightly different 
version of the initial problem, defined by the partition 
function: 



Z sh [V] = Z[V] + Z[-V] 



E 



(c-«') + e W) 



(85) 



which corresponds to a particle in a random potential 
which can explore both V(r) and — V(r). A physical 
realization would be a particle with an Ising spin in a 
random field. As it turns out the physics of this disor- 
dered model is very similar to the original problem. At 
low temperature, it is now related to the distribution of 
the minimum of — |V(r)|. 

We define the generating function of this model 
G s h{x) = (exp(— /j,Z s h[V])), with /i = e l3x , which is re- 
lated to the distribution of the free energy of the particle. 
In the continuum limit and in d = 2, it can be rewritten 
as: 



G sh {x) = H sh [fi] = J DVe- 3 ^ 
S sh [V] = 



+oo /*+oo 



1 



dxd y ( — (W) 2 



(86) 
2 j ucosh(/?y(r)) 



i.e the partition function of the sinh-Gordon model in 
d = 2. Similarly, the d — 1 version of our model is re- 
lated to the well studied boundary sinh-Gordon model 
H defined as: 

G shB {x) = H shB [tj] = J DVe~ s '^ (87) 
SshB [V] = J^°° dy J* dx f^L-(W) 2 

+2^cosh(/3U(a;,0))^ (88) 

Indeed one has, as required, that (V(x, 0) — V(x', 0)) 2 ~ 

4ct In at large \x — x'\, and one only studies (bound- 

ary) observables defined at y = 0. 

In the limit of (3 = +oo one has in both case: 

G sh {x) = Proba(x < mm{V r , -V r )) (89) 
= Proba(x < -max,.|l/,.|) (90) 

and thus the (properly discretized) partition function of 
the (boundary) sinh-Gordon model becomes related, in 
that limit, to the distribution function of the maximum of 
the set of positive random variables |V(r)|. The results 
described in the previous Sections about the statistics 
of extrema of such variables imply that some transition 
must occur as a function of (3 corresponding to a related 
"change of behaviour" in the sinh-Gordon and boundary 
sinh-Gordon models as well. This is a prediction, as we 
are not aware of such a change of behaviour at (3 = /3 C be- 
ing mentionned in the literature. As we now discuss, ex- 
amination of known results is perfectly compatible with 
the transition at (3 = f3 c . 

Let us first describe the known exact results both in 
d = 2 and d = 1. The extensive free energy of the bulk 
sinh-Gordon model is defined as: 



fs 



lim - 

L — >-\-oo 



-L- 2 \nG sh 



(91) 



where the model defined in ( p7| ) is considered in finite 
size L. The model is studied usually using the field 
4> = V\/2/a, the nonlinear term being 2fi cosh((3V) = 
2/j, cosh(bcj)) and its free energy depends on the single vari- 
able b = f3\J a j2 = (3//3 c , where f3 c = \fdj~o is dimension 
dependent. Using the variable b, its exact expression, 
proposed in Ref. [Esl , reads when explicitcd |6ll : 



fM = c 2 {b)^ 

C 2 (b) - 



(92a) 



2tt 



r[ 



2+2b 2 



7rr[l + b 2 ] 

-r[-& 2 ] 
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These results are a priori only valid for b < 1 (\b\ < 1), as 
they were obtained in |5^] from an analytical continua- 
tion of the sine-Gordon model (performing /i — > — // and 
b 2 — ► — 6 2 , A/ being the soliton mass). The constant \i 
was defined in the continuum model by fixing the normal- 
ization of the field < cos(60(r)) cos(6</>r') >= \\x — y\~ 4b 
of the sine-Gordon model. 

The d = 1 version corresponds to the boundary sinh- 
Gordon model usually studied using the cf) = V/^fo 
and 2/it cosh(/?y ) = 2fi cosh(6</>), with again b = (3/(3 c 
[fi c = The analogous expression for the free en- 

ergy reads, from |M: 



fahB(fJ-) = lim -i InGsfes = Ci(&)^ 

L — ;-+og 



(93) 



Ci(b) 



1 + 26 2 



2b 2 



2b 2 



2n 



T[-b 2 } 



Let us now comment on these results. The power law 

dependence in fj, of the free energy is just the naive dimen- 

i 

sional result ~ in both cases. This result should 

hold for (3 < (3 C . However, there is clearly, in both d = 2 
and d = 1 cases, a singularity as (3 — > (3~ as the ampli- 
tude C(b) diverges as b = (3/j3 c ^ 1 ■ This is thus in 
perfect agreement with the existence of a phase transition 
in the particle model. In the sinh-Gordon model itself, 
we do not expect strictly speaking a phase transition, as 
the model is massive both below and above 6 = 1, how- 
ever we do expect some "change of behaviour", which 
may be related to a change of nature of the excitations 
around the ground state. This is not ruled out by ex- 
act results [p2| as it clearly comes here from the physical 
mass acquiring a nontrivial dependence in the bare mass 
parameter /i (contrarily to sine-Gordon model, for sinh- 
Gordon model there is no presently known exact solution 
of a lattice version). 

Let us now interpret these results for our model. They 
mean that the generating function G s h (x) of the free en- 



ergy distribution, with /i 
a traveling wave: 



e^ x , takes indeed the form of 



G sh ~ exp 



with I — hxL and a velocity: 



c = - + a 13 

P 



g(x + cl + 7) (94) 



(95) 



This is exactly the velocity given by the KPP equation for 
the particle model, in the high temperature phase. It also 
yields a front g(y) = cxp(—e ay ) with a = jTprgTg p an d 

7 = i^Tpm p \~nCd(/3/f3c)- This form however should be 
taken with caution as strictly speaking formula (Q) is 
valid only in the limit where L goes to infinity first (at 
fixed fi = e@ x ). It should be compared with the asymp- 
totic behaviour of the front in the region of large positive 
y. We expect universality in the other region of the front 



(of very negative y i.e. x « cl) and exact knowledge 
about this region would be equivalent to exact knowledge 
of the sinh-Gordon model at finite size, which is not yet 
available. 

The physics of the problem of the particle in the ran- 
dom potential leads us to conjecture that the 2d sinh- 
Gordon model (as well as the boundary sinh-Gordon 
model) will exhibit a change of behaviour, the algebraic 
/i-dependance of its free energy will freeze for (3 > (3 C , 
which corresponds to the low temperature glassy phase 
of the particle model. We thus expect: 



1 



(96a) 

(3 < (3 C = y/d/a (96b) 
(96c) 




and presumably log corrections (at least at (3 = (3 C , and 
maybe for all (3 > (3 C ). 

This is confirmed by a renormalization group analysis 
directly on the Sinh-Gordon and Liouville models dis- 
cussed below. 



B. Relation with the Liouville model in d — 2 

The relation between our original model (Q) of the par- 
ticle in the random potential and the Liouville model 
proceeds via the generating function: 



G{x) = (ex V (-e^Z[Vj)). 



exp(- 



V(r)h 



(97) 



which encodes the full probability distribution of the free 
energy of the particle. In the case of the d = 2 potential 
with logarithmic correlations it is identical to the parti- 
tion function of a Liouville model, which one can write 
either on the original lattice or in the continuum (with 
UV and IR cutoffs a and L) as (fi — e@ x ) : 



G{x) = H\p] = / DVe- s[v] 



S l V ]= I d 2 r[—(VV(r)y+f,e 



(98) 



where the functional integral is normalized such that 
H[fi = 0] = 1 (equivalently one redefines H[pi\ — > 
H[fi]/H[n = 0]). We call it the Liouville Model (LM) 
to distinguish it from the continuum Liouville field the- 
ory LFT whose definition is recalled below. A relation 
also exists between the correlation functions of the Gibbs 
measure and some correlation functions in the Liouville 
model: 



<p(ri)..p(r n ) >= 



-/3(y(v 1 ) + ...V{v n ))-S[V] 



(99) 



M>0 
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Strictly speaking the model fl98| ) above is not well defined 
because of the zero mode V(r) — > V (r) + w and must be 
complemented with boundary conditions. In the particle 
problem studied here we have chosen periodic boundary 
conditions with the additional constraint ^ r V(r) = to 
pin the zero mode. 

On the other hand, many results are known for the (re- 
lated) continuum Liouville field theory (LFT), of great 
interest in quantum gravity p3|-|66|,?]. It is usually de- 
fined on a arbitrary genus h manifold with background 
metric g and associated curvature R by the action |67[ : 



S. 



LFT 



1,1 ' ^ 9 ^f 



2b<f> 



47T 



Ry/g</> (100) 



in conventional notations, with (j) = —Vj\f2a and b = 
f3//3 c . The (standard) choice is Q = b+l/b for which the 
theory is critical and has local conformal invariance (with 
a central charge cl = 1 + 6Q 2 = 25 + c) p7[ . It can also 
be formulated as the theory of (liquid) random surfaces 
p6| , |68| e.g. as random triangulations. There one defines 
the total area A = J r e~ 2b, ^ r ) which is nothing but the 
partition function A = Z[V] of the particle problem, and 



which is 



studies the distribution Z(A) 
nothing but P(Z). 

The particle model allows us to make precise state- 
ments on the Liouville model defined above. The LFT 
allows for exact calculations (e.g. of correlation func- 
tions) and in principle one could hope to translate those 
in the particule model (Gibbs measure correlations). The 
relation between the two however is rather subtle. For 
instance, the boundary conditions chosen in the particle 
problem would correspond to Liouville on a torus h = 0, 
except that the additional pinning condition spoils it. We 
will thus not explore here all these intricacies but give a 
few general remarks, mostly about the behaviour of the 
Liouville model under coarse graining. 

First we know that Gi(x) satisfies a RG equation of 
the KPP type. Thus upon coarse graining (i.e as a func- 
tion of I = ln(L/a) see Section III) the Liouville model 
partition function satisfies a KPP non linear RG equa- 
tion. The corresponding front velocity gives the scaling 
of the partition function with [i. The glass transition, 
with freezing of the front velocity, corresponds exactly in 
the Liouville model to the transition between two regimes 
(the so-called c = 1 barrier in the LFT) : 

(i) weak coupling Liouville: b — f3/(3 c < 1. In 
that regime there is no problem to define a continuum 
limit. The KPP non linear RG front solution of velocity 
c(p) = 2(i + js) = c{b) yields Gi(p) ~ ^/U+b 2 ) and in 
the Liouville model the scaling dimensions are the one 
obtained by naive dimensional counting, e.g. : 



d 2 re -PV(r) 



fie 



\nZ[V] 



2{l+b 2 ) 



(101) 



In the LM only one of the two branches (deduced by 
b — > 1/6) is selected as I — > +oo. 



The regime b < 1 is the one where the continuum LFT 
is well defined: there the role of the Q term in (10C) 
is to shift the conformal dimension of the fields e 2 * to 
A(a) = a{Q — a) (while the naive power counting di- 
mensions in the LM are A (a) = — a 2 , see above) which 
renders the Liouville term J e 2h ^ exactly marginal (and 
thus the theory critical). It was argued in |Q that the 
LFT gives the correlations of the Gibbs measure, the op- 
erator p(r) (eq. (g)) corresponding to the Liouville field 
e 2h ^ (and is thus of conformal dimension A (b) = 1 (i.e 
p(r) ~ L~ 2 ). A hint in favor of this conjecture was that 
the corresponding LFT conformal dimensions of the com- 
posite fields p(r) 9 ~ e 2qh ^ is simply A(q) = q(l + b 2 -b 2 q) 
which correctly reproduces the multifractal spectrum: 



J d 2 r p{v) q ~ L 2(1 - 



(102) 



given in Jll| and Section VI (in the weak disorder regime 
q < q c ). This is not a very strong test since the same 
multifractal spectrum can also be obtained within the 
LM model by considering the dimension of the normal- 
ized Gibbs measure (rather than the unnormalized one 



-pv 



). Indeed the effect of the Q term is to shift: 



£ 2bq(p _^ j2bqQ e 2bq<j> ^ Z[V] 



~q e 2bqrt> 



(103) 



To convincingly establish the conjecture of [Q the effect 
of the additional Q term should be checked on the many 
point correlations J67{ , where it is rather more subtle, and 
further investigation is needed. In particular, the RG de- 
scribed here suggests by extension that the Gibbs mea- 
sure correlations (or at least some limits of them) should 
also be computable within the DPCT model. This sug- 
gests a direct relation between LFT and DPCT, a check 
of which would be of great interest. Note also that a 
critical model, which mimics the effect of the Q term (as 
adding an average value to <f> [p7| ) is studied in Appendix 

0.. 

(ii) strong coupling Liouville: b — /3//3 c > 1- This 
corresponds to the glass phase for the particle which, in- 
terestingly, has a non trivial structure. As is well known 
there are serious difficulties in defining the continuum 
LFT in that regime. Using what we know from the par- 
ticle problem we can gain some idea of what happens in 
the Liouville theory. First let us note that since G(x) 
is such that for (3 = +oo and fixed L it is equal to the 
distribution function of the minimum of the set of V(r): 



G(x) = Proba(x < minU(r)) 



(104) 



The (infinitely strong coupling) Liouville model can be 
recast as an extremal statistics problem in that limit. 
The partition sum Z[V] of the particle model being 
dominated, for b > 1 by a few regions of space where 
V(r) ~ Vrain (with little dependence in /?, i.e the Liou- 
ville wall become a hard wall for all (3 < (3 C with thick- 
ness of order O(l)) we expect this spatial heterogeneity 
to show up in LM as well. From what we have learned 
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in previous Sections we know that upon coarse graining 
the following effective Liouville model action S e t / is gen- 
erated: 



Gi(x) = HM = / DV < e~ s ^^ > P 



(105) 



S ef f,l[V} = y d 2 r(^ — (Vf(r)) 2 + z(r)e-^WJ (106) 

i.e a new field z(r) is dynamically generated, and has 
short range correlations but has a broad power law dis- 
tribution: 



that the short scale instability in Liouville may be re- 
lated to the generation of strong inhomogeneities in the 
Gibbs measure p(r), analogous to structures discussed 
in the context of replica symetry breaking. Thus, if the 
mapping onto the LFT is confirmed it suggests to also 
investigate RSB type effects in strong coupling LFT. 



C. Direct renormalisation group analysis of 
sinh-Gordon and Liouville models and traveling 

waves 



P{z)dz ~ z~ 1+ t> 



(107) 



while V = V y is the smooth field introduced in ([L7|) . For 
b < 1 this dynamically generated local field can be aver- 
aged out without changing significantly the action (note 
that even for b < 1 it changes properties of operators 
e -iP v for q > q c ) while for b > 1 it changes crucially the 
physics. One can define the effective Liouville potential 
U[V] for the smooth field V after averaging over the z 
field as: 

U^V] = - In < exp(-ze-^) > Pl{z) 

the bare Liouville potential being U[V] = fi exp(— j3V). 
We can now use the front solution of the KPP equation 
(i.e the scaling region in the large L/a limit) described 
in previous Sections. For b < 1, since < z >p( 2 )< +oo 
we have that for large V 



Let us now illustrate how one can see explicitly the 
freezing of the free energy exponent in the strong cou- 
pling phase from renormalisation group approaches di- 
rectly on the sinh-Gordon and Liouville models. Such 
functional RG methods have been applied to study the 
analogous problem |2j| of the wetting of an interface of 
height V. Related exact RG methods, with various trun- 
cation schemes, have also been applied to the Liouville 
model, and in the context of quantum gravity to the LFT 
|u9f . In all cases we will illustrate how the main physics 
lies in the selection mechanism for the traveling wave so- 
= — hxGt(x = —V) (108^ions of the non linear RG equation. 



U t [V] «c^e 2 ( 1+b2 >*exp(-^0 



(109) 



and thus the coarse grained potential is similar to the 
bare one. However, for b > 1 one has for large V 



Ui[V] « c/j,e 4l Vexp(-/3V) 



(110) 



because of the broad distribution of the z field. 

Since the z(r) are highly heterogeneous on short scales 
it is not surprising that a continuum limit is hard to ob- 
tain for b > 1. These heterogeneities are linked to the 
structure of the glass phase reminiscent of replica sym- 
metry breaking. It is tempting to conjecture that it may 
also be related to the branched polymer structure which 
appear in LFT for b > 1, i.e beyond the c = 1 barrier 
p6[ , or to the spike instability |^58j of fluid membranes. 

Furthermore, let us notice that the LFT theory at 
b = 1 is known to have two marginal operators whose 
dimensions are degenerate e 2f "^ and 0e 2b< ^. This is in ex- 
act parallel with the behaviour of the KPP front solution, 
which develops at b — 1 two degenerate linear eigenmodes 
exp(-/3V) and Vexp{-0V). 

Thus we have seen that the Coulomb gas RG can be 
used to understand the behaviour of the Liouville model. 
A scenario is obtained where for b > 1 new short scale 
degree of freedom are generated (short scale instability). 
Averaging over these changes the effective Liouville po- 
tential. The parallel with the particle model suggests 
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FIG. 13. Liouville wall moving under RG as a trav- 
elling wave. Represented is U[V] on the negative V side, 
of original form U[V] = e _,3 ' v_:c ' ) and also (dashed line) 
G(x) — exp(—U[V]). Both move under RG forming a trav- 
elling front, whose velocity determines the "free energy" ex- 
ponent. For sinh-Gordon a second, mirror image wall is also 
moving symmetrically towards 0. Freezing in the front veloc- 
ity occurs at and below the transition at j3 = /3 C . 

The study proceeds as follows. We consider 

G{x) = H\pL = e? x ] = J DVe~ s[v] (111) 
with S[V] = J d 2 r (^-L-(X7V(r)) 2 + U[V]j 

One can perform a Wilson RG analysis (or if one prefers 
a suitably truncated exact RG analysis), and one finds 
in d = 2 the flow for the local part Ui [V] as: 
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diU = 2U + oU" + 0(U 2 



(112) 



There may also be corrections to a to 0(U 2 ) (in the Sinh- 
Gordon model) , but we focus for now on the RG to linear 
order. Let us recall that the initial condition is Ui = q[V] = 
fie-P v for Liouville and Ui =0 [V] = fie~P v + ^ v for 
sinh-Gordon, and that we are interested in the small fi 
limit. In this limit the initial condition corresponds to a 
a very wide well U[V] (e.g. in the sinh-Gordon model) 
with a very small curvature U" [0] . To obtain the free en- 
ergy exponent as fi — > 0, one simply iterates the RG until 
U"[0] ~ 0(1) at a scale I* — ln(L* /a ) (more precisely 
U"[Q] ~ l/(cao) where ao is the bare UV cutoff of the 
model). At this scale, the free energy is 0(1), as can be 
estimated from gaussian fluctuations (straightforwardly 
at least in the SG model) and thus the initial free energy 
is: 



As we now show, ( |115| ) is incorrect for (3 > (3 C . This 
is so for a subtle reason, as apparently the statement 
that exp( — @V ) is an exact eigenvector of the linear RG 
(and of ( 114 )) cannot fail ! However, by now we are 
well used to fronts: in fact we have encountered ex- 
actly the same equation in our previous solution of the 
REM model via RG (hi(x) in @ is identical to U t {V) 
in (114)). To describe correctly the bare Liouville (or 
equivalently the Sinh-Gordon) model one should gener- 
alize the initial condition Ui=o[V], still assuming that 
Ui = q[V] ~ exp(— (3(V — x)) for V 3> x (x here is very 
negative corresponding to a small /i) . Then one can use 



the saddle point method to estimate (114) as was done in 
( |65| ) to evaluate hi(x) and one discovers that for (3 > (3 C 
the velocity freezes into: 



F 



(113) 



c = 2v / 2ct 
which yields a free energy 



Remarkably, it is now possible to use what we learned 
in the previous Sections and demonstrate the "freezing" 
transition at (3 = (3 C (corresponding to the glass transi- 
tion for the particle) simply from the RG to this order. 
Indeed the solution of the truncated equation is: 



F 



(119) 



(120) 



Ui[V]=e 2 



= [ dV'exp 
iral J 



(V-V'f 



A straigthforward conclusion would then be that the ex- 
act solution corresponding to Liouville is: 



U l [V]=neV + °^ l e-t 3V 
and similarly for the sinh-Gordon: 



since exp(±(3V) are exact eigen vectors of the linear RG 
equation for any (3. From (113) this immediately yields 
the "naive dimensional" result for the free energy 



instead of the naive dimensional estimate, thus inagree- 
ment with our expectation for the SG model ( |96c|) . In 
addition we find that the decay of the renormalized po- 
tential Ui[V] ~ e~ aV is frozen at a — (3 C for all (3 > C 
Ui=o[V\ (114) consistent with ( |110| ). 

What has happened is that although C// = o[V r ] = 
exp(— [3(V — x)) is indeed formally an exact eigenvector, 
it is dynamically unstable, i.e if one chooses another func- 
tion with the same large positive V—x behaviour one gets 
a different velocity (which is not the case for (3 < (3 C ). It 
is easy to see that the choice £/; = o[V] = exp(— (3{V — x)) 
exactly for all V does not make sense for V — > — oo. In- 
deed it is immediately spoiled by the slightest amount of 
coarse graining (as would appear also by considering the 
non linearities in the RG equation) . The simplest way to 
see it is to notice that the coarse grained potential: 



(115) 



(116) 



(117) 



U[V] 



In dv exp( 



-[ie 



-p{V+v) 



V 

2^ 



(121) 



with (3 C = y/2/a. As we know from the above exact 
result this is correct for (3 < (3 C . Note how the poten- 
tial Ui[V] evolves. Using the notation fj, = e l3x (natural 
from our extremal statistics interpretation) it forms a 
"Liouville wall" , which can be seen as a "front solution" 
moving as exp(— [3{V — x — cl)) towards U — (and in 
the sinh-Gordon model there are two symmetric walls 
moving towards U = and reaching it at I = I*). The 
Liouville front velocity is: 



does not grow as ~ exp(— (3V) for large negative V but 
much slower as ~ V 2 . To illustrate further the point let 
us consider the initial condition: 



Ui= [V] 



-0{V-x) 



(122) 



a/3 



(118) 



It behaves as e - ' 3 ^ -1 -' for large positive V — x (and thus 
corresponds to the Liouville model) but goes to 1 on the 
other side. For (3 = +oo it is easy to compute U"[V = 0] 
from (pi since Ui =0 [V] = 6{x - V). One finds: 



which plotted as a function of (3 is the famous parabola, 
such that two values of (3 corresponds to the same c, 
which is also a well known property of Liouville theory. 



U"[V = 0] - e 2l ~^ (123) 



and thus one has that I* defined above is such that: 
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cT 



2V2a 



(124) 



This is in fact valid for all (3 > j3 c as was shown in detail 
in previous sections. 

Thus the freezing transition can be obtained from the 
linearized (i.e lowest order) RG equations, using only el- 
ementary insight from coarse graining or the existence 
of higher order non linear terms. It provides an inter- 
esting example where the naive dimensions hold in some 
regime but are modifi ed in another. Of course, as we have 
seen in Section HID from the study of fronts, to really 
establish the transition and determine the universality 
class one needs to consider higher order non linearities 
in (112) which goes beyond this paper. For the LFT in 
quantum gravity the reader can find some exact func- 
tional RG studies in Ref. p3|. Although not discussed 
in this reference, the non linear RG there seems to also 
exhibit traveling front solutions, whose physics may be 
important in understanding the problem of the c = 1 
barrier. 



VI. CRITICAL DIRAC FERMIONS IN A 
RANDOM GAUGE FIELD 

In this section we relate our RG study of the previous 
section to the study of the critical wave functions of 2D 
Dirac fermions in a random magnetic field. We first con- 
firm the results of fl3f| for the multifractal spectrum, and 
obtain their finite size corrections. Then we study the 
transition from the weak disorder to the strong disorder 
phase, related to the glass transition for the particle, and 
find that the strong disorder phase has a new and non 
trivial structure, leading to what we call quasi-localized 
eigenstates. 



B(r) = 0. The type of correlations studied here cor- 
respond to the most interesting case where the gauge 
potential has short range correlations. In the Coulomb 
gauge, we can introduce the scalar potential <fi such that 
= e^dvej), B(r) = — 9^(r). The gaussian distribu- 
tion of </>(r) is thus given by 



P[<j)] = cte 



x e 



(126) 



where g parametrises the strength of the random mag- 
netic field B. The correlator of the function </>(r) is thus: 



(<Kr)-0(r')) 2 ~2. 9 ln- 



(127) 



In this model, the wave functions at energy E are local- 
ized for all energies other than the critical energy E = 0. 
We restrict our study to the E — critical eigenstate, 
which satisfy: 



HV Q (r) = 



(128) 



For a system of finite size L with appropriate boundary 
conditions there are two independent normalized solu- 
tions of (12S): the first one can be written ^o,i( r ) = 
(#o(r),0) with: 



*o(r) 2 



-20(r) 



£ r e- 20(r,) 



(129) 



the second one being "Jo, 2(1") = (0, VPo( r )) where ^ / o( r )) 
is given by (129) changing <fi(r) — > — 4>(r). We denote ^ r 
having in mind cither a discrete problem, or a continuous 
problem with some short scale cutoff a. 



A. Critical wave function of 2D random Dirac 

Let us first recall the problem of a massless two dimen- 
sional Dirac fermion in a static random magnetic field 
[ p~T] , p~3|j32) . This model, and its non abelian generaliza- 
tions, has received a lot of attention in connection with 
the integer quantum Hall effect transitions with disorder. 
As discussed in p3 two dimensional Dirac fermions can 
experience three generic types of disorder: random gauge, 
random mass and random potential. Random gauge dis- 
order is believed to be a line of fixed points in this gen- 
eral model and is still not yet fully understood. Here we 
address only the random gauge disorder model of hamil- 
tonian: 



H = (ivpdf, - Ap(r)) 



(125) 



where the 01,2 are the 2x2 Pauli matrices and fi = 1,2 
(we set the Fermi velocity Dp = 1 from now on). The 
random magnetic field B corresponding to the gauge 
potential A is chosen to be gaussian with mean value 



B. participation ratios and multifractal spectrum 

Thus in a given configuration of disorder (f>(r) the quan- 
tum probability |\f'o( r )| 2 is identical to the Gibbs prob- 
ability p(r) defined in (|J) for the particle in the loga- 
rithmically correlated random potential V(r) with the 
correspondence: 



l*o(r) 
2<Xr) 



: = p(r) 
(3V(r) 



(130) 
(131) 



and thus the model depends on a single parameter g = 
^(3 2 a. As we have discussed in previous Sections the par- 
ticle in the logarithmically correlated random potential 
undergoes a transition at (3 C — ^/2/cr at which its Gibbs 
measure changes from being dominated by many sites 
(high T phase) to being dominated to a few sites (low 
T phase). Thus in the quantum problem we expect a 
transition at: 



9 = 9c = 1 



(132) 
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with a weak disorder phase for g < 1 and a strong disor- 
der phase for g > 1. In the weak disorder phase the quan- 
tum probability (and thus observables such as the mean 
squared position fluctuations < r 2 > — < r > 2 ) is delo- 
calized (< .. > means averages over ^ )- In the strong 
disorder phase the quantum probability is more concen- 
trated, but it cannot be called localized in the usual sense 
(of an exponential decay around a single center) and in 
fact both phases have rather peculiar properties. 

Properties of wave functions can be described by the 
inverse participation ratios defined from the normalized 
wave function v l'o( r ) m & system of size L by 



R g (L) = / d 2 r\* (r)\ 2q = / d 2 v (p(r)) q 



(133) 



At a very qualitative level, the nature of the eigen- 
function can be inferred from the scaling behaviour of 
the inverse participation ratio with the system size L : 
for an exponentially localized state R q {L) scales [[70| as 
R q (L) ~ const for all q > 0, while for a plane wave de- 
localized state we get R q (L) ~ £-2(9-1). i n addition 
to the localized and delocalized states, there exist states 
such that r(g) = — In R q (L)/ In L is a non linear func- 
tion of q: they correspond to multifractal wave functions 
whose moments cannot be described by a single length 
as usual but rather by a spectrum of exponents. Here, as 
in |l3| ], we also find intermediate multifractal behaviour. 

To compute the finite size inverse participation ratios 
we can use the information of previous Sections since: 



s q {L) = -\nR q (L) = -lnZ q 



q 111 Z 



(134) 



where we have defined Z g = Z(j3 = yj2gja) where Z((3) 
is the partition function for the particle at inverse tem- 
perature (3. In particular we will be interested in the 
multifractal asymptotic scaling exponent r(q) defined by 



r(q) = lim 



(135) 



These exponents were computed previously in |13[| using 
the REM approximation. Here we use our RG results 
and also obtain finite size corrections. Note that these 
correspond to properties of \& defined above and could 
be changed if other boundary conditions were used. 
From the previous Sections we obtain: 



InZg = 2(l + g)lnL + A fi 



,3<1 (136a) 



lnZ 3 = (V?(41nL--ln(lni)) + A g ,g=l (136b) 

lnZ g = ( y /g(4]nL-^\n{lnL))+A g ,g>l (136c) 

where A ff is a sample dependent variable of order O(l) 
with a g dependent distribution (whose tails we have 
characterized previously). From there we obtain s q (L), 
which have different behaviours in the two phases. 



(i) weak disorder phase 

For g < 1 we find, denoting q c = -y= 



s q {L) = 2(q-l)(l-gq)hxL + A q 



2 _ 1 

s q (L) = — (1 -sgn(g)V?) 2 lnL+ -lnlnL + A q , g 



for \q\ < q c 
l - In In L + A 
for \q\ = q c 



s q {L) = 2g(l - sgn(q)^) 2 \nL + -\q\^\n\nL + A q , g 

for \q\ > q c 

where A q ^ g is a fluctuating part of order O(l). 
(ii) strong disorder phase 
For g > 1 we find: 



s q (L) = -2(qV9- l) 2 ln£- -q^/glnlnL + A q<g 

(138a) 

for \q\ < q c 
for q = q c (138b) 
(138c) 



s q( L ) = ~2 In hi £ + 9 
s q (L)=A q ^ g for q > q c 



s 9 (L) = -2|< 7 |Vs(4mL-~mlnL] + A q , g (138d) 



for q < -q c 



s q {L) = -2 ( 41ni- -lnlni 



A n 



for q= -q c 
(138e) 



where A q ^ g is a fluctuating part of order O(l). 

The corresponding scaling exponents r(q) are thus 
identical to the one found in |l3j and in addition we have 
obtained their finite size corrections. In the weak disor- 
der phase one for q > 0: 

2(< ? -l)(l-4 r ) for q<q c = .fl 



r(q) = < ( \ 2 " 7 va (139) 

[ 2q (l - j- c j for q > q c 

which means a parabolic form with a termination point 



at q = q c as represented in the Fig. 14 
T(q) 




parabolic 
approximation 



FIG. 14. Multifractal spectrum in the weak disorder phase 
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In the strong disorder phase g > 1, i.e when q c < 1, 
the above expression becomes (for q > 0): 



r(9) 



-2 1 







y for 9<g c = 

for 9 > q c 



(140) 



Since the inverse participation ratio does not scale with 
the system size L for each integer q, one could naively 
conclude that it is the sign of a localized state (see how- 
ever below). 

As was discussed in Ju| these results can be translated 
into spectrum for exponent a. If one assumes that p(r) 
is of order L~ a in a number L^°^ sites then the above 
spectrum is recovered if: 



/(«) = 8 



[d + — a) (a — rf_) 



(141) 



with d± = 2(1 ± ^fg) 2 for 5 > 1 and d+ = 8y/g, cL = 
for g > 1. It is easy to see that: 



< (r- < r >) 2k >> L 



»((fc+l)/(a)-a) 



(142) 



showing that the eigenstate is never localized in the usual 
sense (exponential decay around a single center) since the 
exponent is always positive for large enough k. Since 



nm« — >±oo max ,mm 

maximum of the Gibbs measure p mt 
the minimum behave for large L as: 



one obtains that the 



Pmax 
Pmin ' 



,£-2(l-^)(l nL )-iv/s 
L-2(i+v^)(i nL )+IV5 



max r p(r) and 

(143) 
(144) 



in the weak disorder phase. 



C. nature of the strong disorder phase: 
quasi- localization 



Let us now concentrate on the case g > 1. There, we 
know from previous Sections that the Gibbs measure of 
t he p article is concentrated in a few sites. Thus from 
(131) the quantum probability |<I>o(r)| 2 is also concen- 
trated in a few sites, analogous to the RSB picture. This 
is a very peculiar type of eigenstate. Indeed if one com- 
putes the quantum average < r 2 > — < r 2 > in a given 
sample, it has a finite probability to be of order 0(L 2 ). 
Thus the eigenstate cannot be considered as localized in 
the usual sense. Since it is peaked around a few sites we 
call it "quasi-localized" . Around these centers the wave- 
function decays fast enough to be normalizable. It would 
be interesting to investigate further the typical spatial de- 
cay of such eigenstates around their (multiple) centers, 
which we expect to be slower than exponential. 



VII. CONCLUSION 



In this paper we have studied the equilibrium prob- 
lem of a particle in a random potential with logarithmic 
correlations, through exact bounds, numerical simula- 
tion, qualitative arguments and a renormalization group 
method that we have developed specifically for this prob- 
lem. We have shown that it exhibits a glass transition 
at finite temperature T c = 1/(3 C > in any dimension. 
This confirms earlier conjectures and allows for a more 
detailed study of the problem. The RG method allowed 
to obtain the universal features of the free energy distri- 
bution at low temperature. The relation to the problem 
of extremal statistic of correlated variables was investi- 
gated. It has been found that it exhibits universal finite 
size corrections, consistent with our numerical calcula- 
tions. 

Most interestingly, we found that this logarithmic 
model provides a particularly simple example (maybe the 
simplest) of a finite dimensional model - i.e with transla- 
tionally invariant disorder correlations- such that the low 
temperature phase is non trivial. It is non trivial in the 
sense that in the thermodynamic limit L — ► +00, there 
are, with a finite probability, several low lying states (i.e 
possible positions of the particle) with energy differences 
of order one, and separated in space by distances of order 
L. Thus the Gibbs measure at low temperature is dom- 
inated by "a few" spatially well separated states. Inter- 
estingly, this transition and this type of glass phase oc- 
curs only for logarithmically growing correlations, faster 
growth (e.g. as in Sinai model) yielding only a glass phase 
with single ground state dominance, while slower growth 
yielding only a high temperature phase. 

Although oversimplified in some respect (it has no in- 
ternal space) it does provide one example of a model 
where the usual droplet picture (which assumes domi- 
nance of a single ground state - or several related by a 
symmetry) does not apply. Rather, it provides one exam- 
ple where some features of the physics usually associated 
to RSB, namely dominance by a few states with exponen- 
tial free energy distributions, can be explicitly exhibited. 
In fact, due to the finite dimensional correlations, there 
are some departures from the behaviour observed in the 
simplest prototype mean field models (such as the REM), 
as can be seen for instance from the free energy distribu- 
tion which has more structure than a simple exponential. 
It would of course be interesting to explore further the 
additional features specific to finite dimensions. 

Although the present model is already of obvious phys- 
ical interest (in 2D it describes e.g. a single vortex in a 
random gauge XY model) its non trivial properties pro- 
vide a motivation to search for models with more degrees 
of freedom and with similar features. One way to proceed 
would be to search for interface models via an internal 
dimensional expansion around the present model. The 
key feature however appears to be the marginality of the 
model, i.e the logarithmic growth of typical energy fluc- 



2G 



tuations. This corresponds to a fluctuation energy expo- 
nent 9 = 0, i.e the situation where the temperature (i.e 
the entropy) is marginal in the RG sense. The droplet 
arguments indeed assume that 9 > 0, consistent with the 
single ground state dominance (and activated behaviour 
typical of a zero temperature fized point where T is for- 
mally irrelevant). In the situation 9 = one does expect 
more generally power laws with T dependent exponents, 
reminiscent of mean field. It would thus be of great inter- 
est to similarly exhibit other non trivial marginal mod- 
els (e.g. spin models with 9 = 0) with similar features 
jnj. Spin models where (domain wall) excitations (in 
root mean square and in average) also scale logarithmi- 
cally (as vortices in the random gauge XY model) are 
presumably good candidates. 

On one hand we have developed a specific (Coulomb 
gas) renormalization group (RG) approach to describe 
the model. From the study of the resulting nonlinear 
(KPP) RG equation, we found explicitly that a freezing 
phenomenon occurs at the glass transition temperature, 
and that in the glass phase a broad (power law) distri- 
bution of fugacities develop - or equivalently exponential 
distribution of local free energy. It is different from more 
conventional perturbative RG (e.g. the one which was 
used to study the dynamics of this model) in the sense 
that the full distribution of probability is followed. This 
turns out to be crucial to describe the low temperature 
phase. 

On the other hand, as we have discussed, two approx- 
imations of the present model, the REM approximation 
and the DPCT hierarchical version can both be solved 
using replica and do require considering the analytical 
continuation to m — ► of contributions of replica sym- 
metry breaking saddle points ]l0| . 

This shows that a RG approach which is explicitly 
replica symmetric but allows to treat broad disorder 
distributions can be consistent with (approximate) ap- 
proaches based on RSB saddle points (j72|. We have il- 
lustrated this on the REM which can be recast in terms 
of non linear RG equations, with a freezing transition. 
In fact one of the striking property of the model is that 
the RG equations derived here are similar - to the or- 
der we have been working - to the one which hold for 
a continuous version of the DPCT problem, the branch- 
ing process. In particular it indicates that both problem 
share the same universal finite size corrections. 

We have also analyzed some connections in 2D (and via 
boundaries in ID) between the model of the particle and 
the Liouville and Sinh Gordon models. The intensive free 
energy of the particle corresponds to the scaling dimen- 
sion in these models with b = (3/f3 c . The glass transition 
corresponds to the weak to strong coupling transition at 
6=1. Beyond, corresponding to the glass phase, the 
scaling dimension freezes as we have also shown via a di- 
rect RG approach on these models. We have seen that 
under coarse graining an additional local field appears in 
the LM and SGM, with broad distribution, and corre- 
sponds to inhomogeneous configurations being generated 



(and broad fluctuations of the local area since the local 
partition function corresponds to local area). 

The present study raises interesting issues to be ex- 
plored concerning the relations with the continuum Li- 
ouville Field Theory (LFT). An outstanding question is 
whether the conjecture of p| is correct for the corre- 
lations. Since we have obtained another result linking 
the problem to the DPCT, the direct comparison of the 
LFT and the DPCT remains to be studied. If it holds 
it means that the conformally invariant many point cor- 
relations can be related in some limits (large separations 
with fixed ratios In r y / In rki ) to the results from the tree 
problem. It would also raise interesting issues about the 
continuation of the LFT beyond 6=1, and its relation 
with the non trivial structure of the glass phase (with 
RSB features) in the equivalent particle model. 

We have also extracted from our approach some con- 
sequences for the problem of the E = critical eigen- 
state of 2D Dirac fermions in a random magnetic field. 
We have confirmed, via our RG method, previous results 
concerning the multifractal spectrum and extracted their 
finite size corrections. We have found that the non trivial 
low T phase of the particle translates into peculiar quasi- 
localized eigentstates for the quantum problem, peaked 
around a few distant centers. It raises the question of 
whether this property can be present in other quantum 
systems. 

Another interesting question is whether the transition 
studied here has a signature in the dynamics as well. 
Note that a similar non trivial structure at low temper- 
ature is also present in the the Sinai model with a bias, 
which renormalizes onto a random walk with algebraic 
waiting times distribution ]5l| . However this is a driven 
system and it would be interesting to see whether non 
driven systems in low dimension can exhibit similar fea- 
tures. 

Finally, an outstanding question is how the present 
model can be studied using 2D conformal field theory 
(CFT) . In particular one wonders what is the signature in 
this context, of the physics which was unveiled here, rem- 
iniscent of RSB, using RG with broad distributions. The 
freezing phenomenon within the non linear RG, which 
transforms the naive scaling dimensions into non trivial 
ones, should correspond to a similar mechanism in CFT. 
Recent progress on CFT classification of disordered mod- 
els where supersymmetry can be used allows to hope that 
such progress is within sight. We hope that the present 
RG method will apply to study other two dimensional 
models with similar features and shed light on the more 
formal field theoretic methods. 

We acknowledge useful discussions with B. Derrida. 
This work was supported in part by NSF grant DMR- 
9528578 (D.C.). 
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APPENDIX A: EXISTENCE OF A TRANSITION 

We use the same method as Derrida and Cook |7^| for 
the directed polymer problem Wm. It is easy to com- 
pute the first two moments pf P\Z], using translational 
invariance and periodic boundary conditions: 



and 



Z = L d e^ r ^ ~ L d +^ CT 

r 

~ BL 2d+2fi ' 2 ° /3<(3 2 



(Al) 

(A2) 
(A3) 



the last estimate being valid as long as the sum over 
r is divergent, i.e (3 < fa = y/d/(2a). The constant 
B > 1 depends on the details of the model, e.g. for d = 1 
one can write B = lirn£,_>+oo f dy exp(— ^-(T^Ly) — 
4<r In L)). Thus for (3 < (3 2 the ratio 



-2 



B 



(A4) 



as L — > +oo. In 
implies that: 



it is shown that the property (A4) 



prob( 



1 



d\nL 



InZ = 



' lnZ)>l/B 



d\nL 



(A5) 



as L — ► +oo. If we take for granted that the fr ee ene rgy 
is self averaging, it implies that for T > T 2 = \j2a jd the 
quenched and annealed (intensive) free energies coincide 
exactly f(T) = f A {T). Thus for T > T 2 the (intensive) 
entropy is s(T) = sa(T) = —drfA(T) and thus one has: 



T>T 2 



(A6) 



Since sa(T) becomes negative below T = T g = y/d/a 
it impl ies that there must be a temperature T c < T 2 at 
which ( |A6| ) breaks down and thus a phase transition. Al- 
though this is harder to prove, it seems that here ( |Arj ) 
holds down to T c = T g . 

Awaiting a rigorous mathematical proof, we have not 
attempted to prove self averaging of /. Not only is it 
highly reasonable in view of our other results but in fact 
if it were not the case, the above argument would imply a 
rather curious - and unphysical - distribution for / (with 
a delta peak of non zero weight smaller than one). In 
addition, as noted in ]7q] , by adjusting the small scale 
details of the model, the constant B can be chosen as 
close to 1 as wanted. 



APPENDIX B: EXTREMAL STATISTICS OF 
CORRELATED VARIABLES 

In this Section we summarize some results on the ex- 
tremal statistics of a set of random variables. We se- 
lected the ones which are useful in putting the problem 



studied here in a broader context. We recall some of 
the classic results from probability theory and we have 
chosen to illustrate them by adding a few simple argu- 
ments which emphasize the importance of some of these 
results to the physics of disordered systems. We denote 
the N random variables either X r , r — 1, ..TV when they 
are normalized in a particular way, or V r when they can 
be readily interpreted as the random potential variables 
studied here (the two differing by a trivial uniform rescal- 
ing V(r) = V r oc X, r ). They apply directly to describe 
d = 1 (N = L) and can be usually extended to d > 1 
(V(r) and N = L d ). 



1. uncorrelated variables 

It is natural to start with the case of N uncorrelated 
variables of identical probability distribution P(V). The 
distribution P(V) can belong to three classes of extremal 
statistics, but we will recall only the Gumbell class. 
Schematically for this class, a well known theorem |25j ] 
states that there exist constants ajv and 6jv such that for 
a fixed y 



Prob{V mm > b N y - a N ) -> cxp(-e y ) 
The constants ajv and 6jv are determined as: 



In 



dVP{V) 



N 



' — oo 

/-a N 
dy I dVP(V) 
-oo J — oo 



(Bl) 

(B2) 
(B3) 



For variables X r chosen from a centered Gaussian of 
unit variance P(X) 
and 6/v as: 



■J=e x / 2 , one can choose etjv 
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b N = 



1 



V21nJV 

a N = V21n/V 1 - ln(47r In N) 

and thus one can then write schematically that: 



(B4) 
(B5) 



X 



min.N 



V21niV- 



l 



V21niV V 2 



ln(47r In N)+y) (B6) 



where y is distributed with the Gumbell distribution 
p{y) = exp(— e v ). 

It is useful to note the property of reparametrization 
associated to a monotonous function ip(V). If one has 
(Bl) for the minimum V m in of the variables V r with the 
constants an d by , one also has (under some weak 
conditions) that (Bl) for the minimum ip(V m i n ) of the 
variables ip(V r ) with the constants a' N — — <xj\r) and 
b' N = bff/t[)'(— ajv)- Note also that we have illustrated 



28 



how to show convergence to Gumbell (and generalized it 
to finite temperature) in the text. 

For completeness we recall the necessary conditions for 
the convergence to Gumbell (i.c P(V) belonging to the 
Gumbell class). First P(V) must decay fast enough at 
V — ► — oo so that there exists yo such that: 



/vo rv 
dy / P(V)dV < +00 
-OO J — OO 

and second, defining 



(B7) 



/?,/,= - ~ L dy L p[v)dv m 



f-oo P (V')dV' J 
one must have for all x < yo 



lim 



t+xR(t) 



P{V)dV 



f-oo P iV')dV> 



(B9) 



These conditions are in fact rather broad. Finally, note 
also the very powerful theorem 2.10.1 of p5| for the rate 
of convergence to the Gumbell fixed point. 



2. correlated variables 



by choosing x = x^, for any a < 1/2. Thus one has 
a general lower bound for the minimum of correlated 
variables. In particular for gaussian variables such that 
V? = 2a In N = 2da In L one gets: 

V mm > -2d^/a\nL + ^/aa\n(Aird\iiL) (B14) 

with probability 1 in the large L limit for any a < 1/2. 
Moreover choosing a = 1/2 and writing: 



V min = -2dyfa\nL+ yfr(- ln(47rdmL) + y) (B15) 



one finds that: 



Prob(y < y) < e y 



(B16) 



This yields a lower bound which can be compared with 
the REM approximation defined in the text. Note that 
the above upper bound is the exact behaviour of the 
Gumbell distribution at large negative y, so in a sense 
the REM approximation saturates the bound in the tails. 
Consequently, to allow for a larger tail (such as ye~ v one 
needs at least a coefficient of In In N strictly larger than 
1/2). 



a. general lower bound 



b. short range correlations and convergence to Gumbell 



We now consider correlated variables with distribution 
P(Vi, ...Vjv). Let us start with a simple but very general 
bound and extract the consequences. One has: 

G{x) = Proba(F mi „ < x) < ^ Proba(F r < x) (BIO) 

r=l,N 

since the reunion of all events V r < x implies the event 
Vmin < x and that Prob(A U B) < Prob(A) + Prob(B) 
(the bound is exactly saturated e.g. when there are 
strong correlations such that V r — V r > > x for all r 7^ r'). 
For variables which have identical one particle distribu- 
tion P(V r ) = fUr^r dV r' p (Vl, —Vn) one has: 



G(x) < N / P(V)dV 



(Bll) 



Let us illustrate the consequences for correlated variables 
X\, ..Xn such that the one particle distribution is a unit 
centered Gaussian. Then it implies for x — > — 00: 



G(x) < 



N 



2ttx 



(B12) 



from which one immediately sees that it implies: 

„ , / r-, ln(47rlnAT) 

Proba X min < x N = - V2 In N + - 



< 



(47rlniV)(2-«) n^+oo 



V21niV 

(B13) 



Let us now consider N centered gaussian variables X r 
with a fixed correlation matrix T rr i = X r X r i, norm alized 
so that T rr — 1. A powerful bound, which refines ( B1C ) 



above, allows to easily demonstrate convergence to the 
Gumbell distribution for a large class of "short enough 
range" correlations. It compares two arbitrary correla- 
tors r£? 7 and r$ with T$ = tH ] = 1. Their associated 
G(x) functions satisfy |25| : 



\G!(x) - G 2 (x)\ < 



.(i) 



->(2)| 



with m rr > = max(\T^,\, |r^,)|. It is obtained by bound- 
ing dG(x)/dr rr > and integrating between I\ and r 2 . It 

will be used to compare T^J, = T rr i with the uncorrelated 

case T rr , = 8 rr ' . 

To adress the question of the universality of the Gum- 
bell distribution, let us now consider a (d — 1) trans- 
lationally invariant correlator T rr r — T(r — r') with 
T(O) = 1 where T(r) is a N- independent function which 
decays to zero as r — r' — > +00. 

Inserting x = a^y — b^ of ( ]B1| ) and (B5) into (B17) 
one easily gets that if T(r) decreases fast enough, one 
has G(x = cbNy — &a0 = exp(— e^) at large N, i.e one has 
convergence to the Gumbell distribution with exactly the 
same coe ffici ents and b^ as in the uncorrelated case, 
so that (B6) still holds. As one sees by studying the 



U/2 



(B17) 
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bound this result holds as long as T(r) decreases faster 
than l/ln(r) (this is theorem 3.8.2. of pft). The lim- 
iting case (which does not satisfy Gumbell, as discussed 
below) being T(r) ~ r/ ln(r) at large r. 

Let us give a simple self consistency argument, more 
enlightening than the bounds, which explains why T(r) ~ 
1 / In r should be the limiting case between the short 
range (Gumbell) universality class and other behaviours. 
Let us split a set of 2 N correlated variables X\ , ..X2N into 
subsystem 1, Xi,..Xn, and subsystem 2, Xn+%, --X^n- 
If correlations are very short ranged (e.g. exponentially 
decaying) it seems reasonable to first neglect correlations 
between 1 and 2 and find the minimum in each subsys- 
tem, which read respectively: 



Y ~ frnCr i HLrhiiV) Xj 

Amin,, « V21niV- + 

2 V21niV V21niV 



(B18) 



with i = 1,2 and where ij, X2 are independently dis- 
tributed with the Gumbell distribution. The symbol V 
indicates that the minimum (in each subsystem) is with 
respect to a slightly different distribution than the origi- 
nal one, since all cross correlations between the two differ- 
ent subsystems have been set to zero. The second stage 
is to add the correlations between the two subsystems. 
Typically, the minima 1 and 2 will be a distance ~ N 
apart and thus their original cross-correlation is ~ T(iV), 
and thus, for short range correlations, very small com- 
pared to the fluctuating part x. k / V2 In N. Thus the dis- 
tribution of the mininum X^ 2 ^ of the original 2N vari- 
ables should be given with better and better accuracy at 
large N, as X^^J = min(X mini i, X mirit2 ) ( which is au- 
tomatically satisfied by the approximation ( B18 )). The 
corrections are irrelevant at large scale provided the typ- 
ical root mean cross correlation between the subsystems 
remain smaller than the typical fluctuations of the mini- 
mum in each subsystem, a condition which reads: 



y/T(N) < 1/vTnJV 



(B19) 



which indeed gives correctly the basin of attraction of 
the Gumbell distribution. Furthermore, in the limiting 
case F(r) ~ r/lnr the above argument shows that the 
distribution of the Xi should be changed, which is also 
the case, as we now examine. 

So, to summarize, if correlations are short ranged with 
r(r) decreasing faster than l/ln(r) this is the "SR uni- 
versality class" . It includes the REM, and one can check 
that the finite size corrections in |(| are reproduced (at 
T = 0). 



V' = v r + u 



(B20) 



with V r a set of uncorrelated gaussian variables, U a gaus- 
sian variable uncorrelated with the V r . The n c learly, if 
one chooses the variance of U big enough (B6) cannot 
hold. To keep using normalized variables (r rr = 1) one 
defines: 



XL — 



1 



^1 + w N 



,X r 



WN 



1 + W N 



(B21) 



where u is a centered gaussian random variable with 
unit variance. The correlation matrix is then r' ; = 
i+ Wlf ($rr' +w N ). Clearly one has: 



X' 



1 



x„ 



W N 



1 + w N 



(B22) 



Using the expression ( |B6| ) for X m i n one sees that for 
deviations from Gumbell to arise one needs that wn ~ 



r/laN. In that case one gets from (B6) that 



X' 



r—— lln(47rln7V) 

V21niV \ ' 

2 V2hi7V 



2tu- 



V21n7V 



(B23) 



These simple considerations thus allow to understand 
simply the li mitin g case, that if T(r) decreases as r/ ln(r) 
one has that (B6) still hold (with the same constants) but 
the distribution of y— t now converges instead to the con- 
volution of the Gumbell distribution and the gaussian of 
variance 2r (see e.g. theorem 3.8.2. of p5|). 

Increasing the range of correlations even further, one 
gets into a regime where the fluctuating part (in the 
X variables) is larger than 1/VhiiV (and thus in the 
V ~ XylriiV variables the dominant finite size correc- 
tions are non sclfavcraging). The fluctuations become 
then entirely gaussian, being controlled by the U part, 
i.e the q = mode. For instance, if T(r) decreases as 
l/(m(r)) a with 1/3 < a < 1 then (theorem 3.8.4. of 
p5|) one has: 



P{V min > -T(N) 1/2 x - (1 - r(A0) 1/2 V2~hTiV (B24) 



(21niV- -ln(47rlniV))) 



r 2TT- 1 / 2 e- x2 / 2 (B25) 
J — 00 



As illustrated below, this behaviour (entirely con- 
trolled by the q = mode) is in a sense more long range, 
and further away from Gumbell than the problem of log- 
correlated variables that we are interested in and that we 
now discuss. 



c. long range correlations and absence of convergence to 
Gumbell 

There is a simple but instructive model of correlated 
variables which can be easily solved and illustrate cases 
where Gumbell does not hold. If one takes: 



d. log-correlated variables 

The case of log-correlated variables is difficult and little 
is known. We just make a few comments. 
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T(r) *lnN 



InN 



(lnNj- a 



Gumbell 



In r 



InN 

FIG. 15. Correlations as a function of lnr. The straight 
line corresponds to the log-correlated variables studied here. 
The thick line corresponds to the limit where the short range 
Gumbell (and REM) behaviour holds, with T(r) ~ l/(lnr) a 
and a < 1, the curved solid line to the case where a convo- 
lution of Gumbell and gaussian holds (marginal case a — 1) 
and the dotted line a > 1 when the mode q — dominates 
the behaviour. 

Let us first discuss the form of the correlator. The cor- 
relator (for the normalized variables X r = V r /y/2a InN 
in d = 1) is of the form: 



r(r) = 1 



lnr 
IdJV 



(B26) 



for 1 ^> r ^S> N. One must then distinguish the two 
other regions. For small r, the precise form could vary 
by adding a short range correlated noise. This is what 
we term short scale details, and an important question is 
the extent of universality of the results (scaling of min- 
ima, distribution) with respect to the small r form. For 
r ~ L the behaviour depends on boundary conditions, 
which may also be important (see below). For the pe- 
riodic system in the simulation T(r) = T(N — r) and 
F(r) actually becomes negative at r = N/2 and of order 
c/ In TV (see Section IV). Adding a small uniform U noise, 
as described above in B2c ( could make T(N/2) = 0, 
so generally speaking one can discuss forms such that 
r(7V/2) = 0. Seen as a scaling functio n of z = -nrj?, r(r) 
then converges for large N towards ( |B26| ), but it does 
have boundary layers at z = and z = 1. 

It is useful to plot on the same graph the various cases 
studied in this Section. This is represented in Figure |l5| . 
We have represented schem atica lly T(r) hiN versus lnr, 
for the log-correlated form ( B26 ) above, and for the vari- 
ous cases r(r) ~ [ETrJ 5 w ^ n a > 1 (Gumbell behaviour), 
a = 1 and a < 1. 

As discussed above in the log-correlated case the be- 
haviour of r(r) near lnr = In AT can be considered as 
uncertain to order 1/lnA. This can be seen either from 
the q — mode, which depending on boundary condi- 
tions one may adjust by this amount, as discussed above, 
or even looking at the first non trivial mode, q = 2tt/L 



which has a contribution of the same order. We know 
from the previous paragraph that these contributions can 
shift the x variable by a gaussian, so it makes it unlikely 
that the Gumbell distribution would hold in that case. 

To conclude, we have given the various behaviours as 
a function of the range of correlations. The presence of 
the In In N corrections seems to be more robust than the 
distribution of x. For the marginal case with q = LR 
disorder, the same \ In In N corrections hold as for the 
REM, while the distribution is changed. On the other 
hand, for log-correlated variables, we expect a different 
coefficient | In In N as discussed in the text (and we do 
not expect the Gumbell distribution to hold). 



APPENDIX C: GUMBELL VIA RG 

In a more detailed analysis Eq. (^30|) yields 
lnln(l/Gz(a;)) = I + lnln(l/Go(x)) which can be rewrit- 
ten in a front-like form: 



Gi(x) = exp(-e l+ ^ x) ) 



(CI) 



where 4>{x) = lnln(l/Go(ir)). In this Appendix we set 
dl — > I, The center of the front is at x — — m(Z) 
solution of (j>{— m(l)) — —I. One can Taylor expand 
4>(x) = —l + y + ^5iy 2 + .. with y = oei(x + m(l)), 
ai = (t>'(-m(l)) and Si = tf' {-m{l)) / (j>' {-m{l)) 2 . Thus 
in the variable y, Gi converges to a Gumbell limit distri- 
bution Q(y) = exjp(—e y ). It holds provided higher terms 
in the Taylor expansion are irrelevant (a necessary, and in 
simplest cases sufficient, condition being that the second 
one Si — * 0). 

If no rescaling of disorder is performed, in the rele- 
vant large negative x region one has Go (a;) ~ 1 and thus 
4>{x) ~ ln(l — Gq{x)). Two cases must be distinguished 
because the limit T — » and N — > +oo do not commute: 

(i) finite fixed temperature T > 0: then for x — > — oo 
one has 1 - G (a;) - GiC^e^l + O(e^)) with C k = 
f v P(V)e~ kl3V and we assume that Cx,C% < +oo exists 
(distributions falling faster than exponentials) . Then the 
situation is simple as (f>(x) = fix + lnGi(/3) + 0{e^ x ) 1 
m{l) - l/P + l/0laCi(0), ai ={3 and 0" (x) / <p' (x) 2 -> 
exponentially fast. For a Gaussian distribution: 



m(l)~ij- 



1 



-a(3 



(C2) 



There is no transition to a glass phase. 

(ii) zero temperature: it is an extremal statistics prob- 
lem. Then clearly 1 — Grj(x) does not decay as an ex- 
ponential. Let us consider a class of distributions such 
that 1 - G (x) - {A\x\)~^ exp(-{B\x\) a ) with a > 1 
(plus exponentially small corrections). This contains the 
Gaussian (of variance a), of most interest here, for a = 2, 
B = l/y/2a, 7 = 1 and A = y / 2Ti/a. Then one easily 
finds from above that 
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mm sb !(/- 2ln/- 7 lnM/ J B))- 
B a 



a, 



Ba(l --\nl~ 1 MA/B)) 1 - 



(C3) 
(C4) 



and that 4>"/<f>' ~ VM" thus the convergence to the 
Gumbell front holds. Note that the quantity aim(l) ~ 
al — j hi I + 0(1) exhibits some universality. 

One thus recovers the standard theorems for extremal 
value statistics reviewed in Appendix [b[ and the relation 
to the normalizing constants defined there as: 



m(l) = aN cti = 7 — I = In iV 
on 



(C5) 



In the Gaussian case, using the val ues given above one 
finds that (C4) indeed yields Eq. B5 in Appendix [§| 
(up to subdominant terms). Although the distribution 
is universal, the normalizing constants obviously depend 
on the details of the tail of the distribution. Note in all 
cases the presence of finite size corrections involving a 
logarithm. 

There is a very small temperature (J3l ~ y/hiL for 
Gaussian) where the behaviour of Gq(x) changes from 
(i) to (ii). It can be seen seen by rescaling temperature 
or equivalently disorder, with system size as in the REM. 

Let us examine the case where the constants Ai and Bi 
are rescaled and now /-dependent (see also e.g. p8[). One 
can still use formulae flC4|). Let us choose Bi = bl~ 1+1 / a 
and Ai/Bi = est (which includes the Gaussian REM). 
One finds at T = that m(l) ~ ±(Z - ^ In I - % hx(A/B)) 
and ai ~ ba. In the gaussian case ai = 2al one recovers 
the REM result: 



m(l) re ^(21 - i 1el(4tt/)) 



(C6) 
(C7) 



at T = (i.e ( f41a| ) setting I -> d/ and a -> cr/d). The 
analysis can be performed at any T and now yields a 
transition temperature when the behaviour of Go,; (a;) at 
large x changes. 



APPENDIX D: AN EXTENDED MODEL 

A richer phase diagram can be obtained by adding a 

Jln^ 



logarithmic background potential |75| V (r 



to the previous random potential (Vrf(r) — Vd(r')) 2 ~ 
Aaln^^- for a < |r - r'| < L and V d {r) = (i.e 
writing V(r) = V (r) + V d (r)) in (Q). The choice of the 
origin breaks translational invariance. The competition 
between the disorder and the binding background poten- 
tial (which if strong enough, tends to favor localizing the 
particle in wells far from r = 0) yields the phase diagram 
of Fig. Another closely related model (model II) 

which preserves statistical translational invariance and 
has the same phase diagram is: 



Z V [V] 



-PJ 



-pV d (r) 



(Dl) 



describes a problem with either zero or one particle 
(vortex) present, the energy cost of the vortex being 
J\n(L/a). It is thus a one vortex toy model of the 
recently studied XY model with random phase shifts 

us 

In the absence of disorder the model with a back- 
ground potential (model I) trivially exhibits a transition 
at P = p* = d/J. At low temperature (5 > (3* the 
Gibbs measure is p(r) ~ C(^) /3J with C = Z L=oc a 
finite constant and the particle is bound to r = (it 
has a finite probability to be within a fixed distance of 
r = 0). At high temperature j3 < (3* the Gibbs measure 
becomes p(r) ~ [j,) d ~^ J {f) l3J arl d the particle is delo- 
calized. This transition can be seen in the free energy 
density f = F/hxL = -T In Z/ In L since: 







-U(3* 



f3>(3* 

0) P<0* 



(D2) 
(D3) 



for (3 < (3* . This first order transition occurs as / reaches 
its bound (since Z > 1 due to the lattice cutoff, one has 
that / < 0). The model II has the same / and a similar 
transition with either one vortex present (3 > (3* or zero 
(3<(3*. 

In the presence of disorder the RG developed in this 
paper can be extended straightforwardly and leads to: 

-AG(x) = -d x G + ~%G + F[G] (D4) 
a ad 

F[G] = -G{1 - G) (D5) 

The additional term thus results in a simple shift in the 
front velocity. The position of the front m(l) thus leads 
to the free energy / = m(l)/(dl) which can have three 
distinct analytical forms: 



(3m(l)/l ~ dpftf) = 

-(d + a(3 2 - J (3) 







high T phase I (D6) 
J (3) localized phase II (D7) 
bound phase III (D8) 



The phase diagram is represented in Fig [l6| using the 
reduced temperature T/J and the dimensionless disor- 
der parameter a = a/ J 2 . For Ada < J 2 one defines 
f3*(a) = ±(J - VJ 2 - Ada). The RG analysis yields 
three phases. In the model with the background poten- 
tial (model I) they are, respectively: 

the high temperature phase (for (3 < /3 C when Ada > J 2 
and for j3 < (3* (a) for Ada < J 2 ): Entropy wins and the 
particle is delocalized over the system. 

the localized phase (for a > a c = l/(Ad) and [3 > f3 c = 
yjd/a): The KPP velocity is frozen. The disorder wins 
and the particle freezes in wells away from the origin. 
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tie bound phase (for a < a c — l/(4<f) and /3 < /?*(g)): 



mat/9906076 



The particule is bound to the origin. Within this phase 
near the phase boundaries (where the bound state length 
is large) a a crossover can be distinguished as a remnant 
of the freezing transition. The bound phase arises be- 
cause of the bound / < (or equivalently the velocity of 
the KPP equation must remain positive). 

In model II the bound phase correspond to no vortex 
present. When one vortex is present it can be either lo- 
calized in a few wells or in a high-T phase (as studied in 
the text of this paper) . 

Both transitions away from the bound phase are first 
order while the transition between high temperature 
phase and localized phase is continuous. An interest- 
ing feature is the multicritical point where the transition 
becomes continuous. 



l/4d 




high temperature 



l/2d 1/d 

FIG. 16. Phase diagram in presence of both disorder and 
external potential. The freezing of the KPP velocity still oc- 
curs at f3 = p c and is represented by the dashed line and its 
solid prolongation: it remains a transition line for a > AdJ 2 
and becomes a crossover line for 4da < J 2 . 
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